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INTRODUCTION 


An  Important  tool  for  analyzing  the  relative  effectiveness  of  gun  air 
defense  systems  is  a plot  of  lines  of  constant  (ISO)  probability  of  kill 
(PK)  in  the  volume  of  space  surrounding  the  planned  deployment  of  the 
system;  hence,  the  name  ISO-PK.  These  plots  are  usually  presented  as 
two-dimensional  sections  taken  at  selected  places  in  the  volume.  The 
lines  of  constant  kill  probability  sometimes  form  eliptlcal  traces  and, 
perhaps,  for  this  reason,  they  are  commonly  referred  to  as  "footprints”. 
This  report  presents  the  computer  program  and  explains  the  algorithm 
which  calculates  the  burst-kill  probability  and  automatically  plots  the 
footprints  on  a computer  line  printer.  Also  Included  in  this  report  is 
an  analysis  of  weapon  pointing  errors  which  includes  the  effects  of  lead 
angle  generation.  This  model,  which  may  be  termed  an  error  budget  model, 
was  originally  developed  for  the  Air  Defense  Evaluation  Board  (ADEB)^. 

The  expression  for  burst-kill  probability  represents  a considerable 
improvement  in  methodology  for  performing  sensitivity  analyses  of  perfor- 
mance parameters  of  certain  gun  air  defense  weapon  systems.  This  expres- 
sion is  combined  into  a computer  program  that  produces  plots  on  a line 
printer,  thus  greatly  reducing  the  amount  of  manual  effort  needed  to 
perform  footprint  studies.  The  program  has  been  used  most  extensively 
for  sensitivity  analyses. 

Sections  of  the  report  describe »in  order,  the  operation  of  the  program 
called  "ISO-PK",  the  error  analysis  which  determines  the  weapon  pointing 
error,  and  the  input  data  requirements  for  the  ISO-PK  program.  The  pro- 
gram logic  diagram,  program  FORTRAN  source  listing,  and  a sample  problem 

^WECOM  ADEB  Task  Force,  Analysis  of  Air  Defense  Gun  Systems.  Technical 
Note  RE-TR-70-191,  Research  and  Engineering  Directorate,  HQ,  US  Army 
VTeapons  Command,  Rock  Island,  IL,  October  1970. 
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are  In  the  final  sections. 

FOOTPRINT  GEOMETRY 

Consider  a coordinate  system  with  an  air  defense  gun  system  at  the 
origin  of  the  axes  as  shown  In  Figure  1.  The  azimuth  angle,  a,  is  measured 
counterclockwise  from  +X  in  the  X-Y  plane.  The  elevation  angle,  e.  Is 
measured  upward  from  the  X-Y  plane  to  the  line  from  the  origin  to  P (this 
line  Is  called  the  slant  range , r) . 

Let  the  X-Y  plane  represent  the  ground  plane  and  z represent  altitude. 

The  target  is  assumed  to  be  at  a point  in  space  with  velocity  components 
(x,  0,  0)  (the  dot  notation  here  refers  to  the  first  time  derivative  of 
the  variable).  It  is  of  particular  Interest  to  know  how  certain  perfor- 
mance characteristics  of  the  gun  affect  its  overall  performance  in  tracking 
and  firing  at  a given  target  at  points  along  straight  lines  parallel  to 
the  X-axis.  These  lines  may  be  at  various  altitudes  or  at  various  values 
of  Y (crossing  ranges). 

The  air  space  around  the  gun  is  divided  into  a 3-dlmenslonal  grid  so 
that  only  a discrete  number  of  points  P need  to  be  considered.  The  user 
may  select  any  plane  parallel  to  the  X-Y  plane  (Z  ^ 0)  or  any  plane 
parallel  to  the  Z-X  plane  (Y  ^ 0) , see  Figures  2 and  3.  Given  the  volume 
of  space  to  be  considered  in  terms  of  limits  on  X,  Y,  Z,  the  ISO-PK  program 
establishes  a grid  in  two  quadrants  of  the  selected  plane.  The  grid  con- 
tains 100  Increments  horizontally  and  50  increments  vertically,  corresponding 
roughly  to  the  size  of  one  page  of  computer  paper  on  a printer  with  spacing 
increments  of  1/10"  by  1/6",  respectively. 

In  the  plane  there  are  51  x 101  ■ 5151  points  denoted  by  Pij ; i=l,2,***, 
101;  j“l,2,‘**,  51.  If  the  plane  under  consideration  is  parallel  to  the 
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Figure  3.  Gun-Target  Geometry  - X-Z  Plane 
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X-Y  plane  (Z  ^ 0),  the  Y position  refers  to  the  crossing  range,  (see  Figure 
2).  If  the  plane  is  parallel  to  the  Z-X  plane  (Y  ^0),  the  Z position 
refers  to  the  altitude  (see  Figure  3) . 

• • • 

The  ISO-PK  program  computes  the  variables  x,  y,  z,  a,  e,  r,  a,  e,  r, 
at  each  point  . The  gun  system  errors  are  determined  as  a function 
of  these  variables  from  gun  system  performance  data  which  is  provided  as 
input.  The  resulting  gun  system  errors  are  used  to  compute  the  burst- 
probability  of  kill  (BPK)  from  the  analytic  formula  in  reference  2.  The 
BPK  is  calculated  about  the  point  of  predicted  intercept  of  the  target  and 
the  projectile.  However,  the  value  of  the  BPK  is  affixed  to  the  point 
corresponding  to  the  present  position  of  the  target  (P^j ) . Certain  values 
of  BPK  (e.g.,  0.05,  0.1,  0.5)  and  corresponding  characters  (e.g.,  H,  I, 

K)  are  specified  in  the  program  to  correspond  to  the  constant  BPK  lines 
printed  for  the  footprints.  A character  corresponding  to  the  numerical 
value  of  the  BPK  is  printed  whenever  the  computed  value  of  BPK  at  a point 
P^^  equals  or  crosses-over  one  of  the  specified  values.  The  lines  of 
constant-kill  probability  then  appear  as  traces  of  printer  characters  in 
the  plane  of  the  paper.  The  matrix  of  computed  BPK  is  saved  during  the 
horizontal  raster  scan  and  is  scanned  again  in  vertical  order  to  fill  in 
the  contour  lines  and  eliminate  any  apparent  discontinuities  in  the  printed 
lines  when  the  lines  have  infinite  slope. 

Two  separate  computations  are  performed  at  each  point  P^^  and  the 
results  of  both  are  superimposed  on  the  computer  plot.  The  first  compu- 
tation performed  is  the  BPK  determination  as  previously  discussed.  The 
2 

Banash,  Robert  C. , An  Analytic  Procedure  for  the  Computation  of  Burst- 
Kill  Probabilities  in  Air  Defense.  Technical  Note  SY-TN9-70,  Systems 
Analysis  Directorate,  HQ,  US  Army  Weapons  Command,  Rock  Island,  IL, 

October  1970. 
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second  is  a determination  of  possible  dead  zones  caused  by  gun  system 
constraints*  The  azimuth  and  elevation  rates  of  the  gun  are  computed 
and  compared  with  maximum  values  allowed  as  specified  in  the  input  data. 

If  either  are  greater  than  the  maximums  allowed,  the  letter  A (azimuth) 
or  E (elevation)  will  be  printed  corresponding  to  the  type  of  constraint 
violated  (if  both  are  violated,  only  an  A will  appear).  The  gun  elevation 
is  also  calculated  and  compared  with  the  maximum  elevation  allowed.  If 
the  limit  is  exceeded,  an  E is  printed.  The  program  also  prints  the  boun- 
dary formed  by  the  maximum  effective  range  of  the  gun.  The  range  which  is 

used  in  this  comparison  is  the  range  at  predicted  intercept  and  not  the 

slant  range  at  time  of  burst.  The  character  R is  printed  to  signify  vio- 
lation of  this  constraint.  Next,  at  each  point  a hypothetical 
situation  concerning  total  system  reaction  time  is  evaluated.  The  letter 
T is  printed  if  the  target  could,  at  its  present  velocity,  move  outside 

the  range  of  the  weapon  if  it  was  first  detected  at  that  point.  This 

provides  a means  of  making  a comparative  evaluation  of  a system's  ability 
to  initiate  an  engagement  against  a target  suddenly  appearing  within  the 
lethal  volume  of  the  system.  In  any  particular  problem,  the  analyst  may 
also  insert  additional  constraint  conditions  into  this  part  of  the  program. 
ERROR  ANALYSIS 

The  BPK  subroutine  used  in  the  ISO-PK  program  requires  the  means  and 
standard  deviations  of  the  errors  associated  with  the  gun  system.  These 
errors  are  divided  into  two  classes.  One  is  the  within-burst  errors. 

These  are  the  errors  associated  with  the  projectile  (e.g.,  the  residual 
dispersion  and  muzzle  velocity  variation) . The  other  class  of  errors  is 
the  burst-to-burst  errors  due  to  inaccurate  gun  pointing. 
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The  gun-pointing  errors  can  be  further  subdivided  into  two  categories. 
One  la  the  error  in  the  gun  position  due  to  errors  in  sensing  the  present 
position  of  the  target.  It  is  assumed  that  a gun  system  will  measure  the 
present  position  of  the  target  and  then,  based  upon  these  measurements, 
determine  the  lead  angle  necessary  to  achieve  a hit  on  the  target,  given 
a certain  open  fire  time.  If  these  measurements  are  in  error,  the  system 
will  produce  an  Incorrect  lead  angle  (a  gun  pointing  error). 

The  other  source  of  gun  pointing  error  is  the  gun  servo  mechanism. 

This  error  is  due  to  the  inaccuracies  which  are  usually  present  in  the 
servo  system  which  positions  the  gun  at  a generated  lead  angle.  Figure  4 
shows  a schematic  of  a gun  air  defense  system  with  the  various  error 
sources  noted  (the  "Ballistic  Errors"  are  the  within-burst  errors).  The 
program  will  accept  the  mean  (bias),  the  standard  deviation,  and  the 
covariances  of  the  random  components  of  each  of  the  error  sources. 

An  analysis  of  the  gun  pointing  errors  was  reported  in  reference  2. 
This  analysis  assumed  that  sensor  errors  produced  gun  pointing  errors  of 
equal  magnitude.  Also,  it  assumed  there  was  no  lead  angle  generation. 
Consequently,  the  time  of  sensing,  the  time  of  fire,  and  the  time  of  pre- 
dicted Intercept  were  all  the  same  (no  time  of  flight).  A modification 
to  this  analysis  is  presented  here  which  includes  the  time  of  flight  and 
the  effects  of  lead  angle  generation  and  approximates  the  lead  angle  error 
due  to  sensor  error.  The  azimuth  and  elevation  errors  produced  from  this 
analysis  is  combined  with  the  gun  servo  errors.  These  combined  errors  are 
then  projected  into  a two  dimensional  plane  at  the  target  at  the  time  of 
predicted  intercept.  Note  that  the  target  position  at  predicted  inter* 
cept  is  not  the  same  as  the  present  position  of  the  target. 

^Ibld. 


10 


Sensing 

Errors 


(r. 


/ 


^Senso^^ 

I 


Lead  Angle 
Servo 
System 


I 

II 

T 

I 

I 

I 


I 

I 


X 


¥ 


Lead  Angle 
Computation 


Gun  Position 

^5erv?T?offlS!?3?7 


Future 

Target 

Position 


Miss 
Distance 


Proj ectile 

Gun  Servo 

Proi . 

Ballistics 

^ Weaoon 

t"  PTtl 

Position 

^Pointing  Position 

oy  oim  cui 

Ballistic  Gun  Servo 

Errors  Error 


Figure  4.  Schematic  of  the  Basic  Components  of  a Gun  Air  Defense  System 


Sensor  Errors 


To  determine  the  gun  pointing  error  due  to  sensor  errors,  a first- 

order  error  analysis  is  performed  on  the  equations  defining  the  position 

of  the  gun.  The  independent  variables  in  this  analysis  are  the  measure- 

• • • 

ments  of  the  present  position  of  the  target  (r,  a,  e,  r,  a,  e) , and  the 
dependent  variables  are  the  azimuth  and  the  elevation  of  the  gun. 

Define  the  gun  barrel  position  by  the  vector  Z,  (Z(l)  = azimuth, 

Z(2)  = elevation),  and  the  vector  of  functions  which  define  Z to  be  f(X). 
This  is 

Z - f(X) 

Where  X is  the  vector  of  measurements  taken  at  the  present  position  of 
the  target 

X = (r,  a,  e,  r,  a,  e)^  . 

The  functions  which  define  the  gun  position  are  derived  in  the  Appendix. 

The  first-order  approximation  of  the  change  in  the  dependent  variables 
(AZ)  resulting  from  a small  change  in  the  Independent  variable  (AX) , 
about  a nominal  point  X = Xo  can  be  expressed  as 

Az  = 3f^X)  AX  . 

The  small  change  AX  may  be  Interpreted  to  be  a random  error  associated 
with  the  measurements  of  the  present  position  of  the  target.  That  is, 

AX  is  a random  variable  with  mean  and  covariance  matrix  COV (X) . 
Therefore,  AZ  is  also  a random  variable  describing  the  errors  in  the 
gun  pointing  position.  Define  the  mean  value  of  AZ  to  be  and  its 
covariance  matrix  to  be  COV(Z). 
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Let 


F.iMl 

3X 


which  Is  the  Jacobian  matrix  evaluated  at  the  nominal  value  , Then 
the  mean  value  of  AZ  is 

E [az]  = F ECAX] 


or 

M2  = F Mx 

£[•]  denotes  expected  value  of  the  random  variable.  Also,  the  covariance 
matrix  of  the  errors  AZ  is  given  by 

COV(Z)  = F COV(X)  f"^  . 


The  partial  derivatives  necessary  to  define  F in  FORTRAN  notation  are 
listed  in  the  section  titled  FORTRAN  Program  Source  Listing,  page  56. 

The  statistics  of  the  measurement  errors  (M^  and  COV(X))  are 
assumed  to  be  constant  (in  time).  However,  provisions  have  been  made  in 

the  ISO-PK  program  for  a user-supplied  subroutine  which  can  modify  these 

• • « 

errors  as  a function  of  the  engagement  parameters  (a,  e,  r,  a,  e,  r) • 

This  subroutine  is  called  MESIG. 

Gun  Servo  Errors 

The  gun  servo  errors  will  be  considered  next.  Define  the  mean  servo 

error  to  be  the  two-vector  M and  the  (2  x 2)  covariance  matrix  for  the 

s 

servo  errors  to  be  COV(S).  The  first  row  in  and  COV(S)  corresponds 
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to  azimuth  while  the  second  row  corresponds  to  elevation  error.  The  mean 
values  of  the  gun  servo  errors  are  assumed  to  be  constant  for  a gun  system. 
These  bias  type  errors  may  be  likened  to  a bore  sight  error.  The  covari- 
ance matrix  Is  assumed  to  be  a function  of  the  gun  rates.  It  Is  also 
assumed  that  the  azimuth  and  elevation  errors  are  Independent  of  one 
another,  so  the  off-dlagonal  terms  of  the  matrix  COV(S)  are  zero.  The 
standard  deviations  of  the  errors  on  the  azimuth  and  elevation  directions 
are  Input  to  the  ISO-PK  program  as  a table  of  values  which  are  a function 
of  gun  rates.  A corresponding  table  of  gun  rates  Is  also  Input  to  the 
program.  In  the  program,  the  appropriate  standard  deviations  are  deter- 
mined by  looking  up  the  rates  In  the  table  and  linearly  Interpolating 
between  the  tabular  values.  They  are  then  squared  to  obtain  the  variances. 
It  Is  expected  that  detailed  system  performance  models  of  the  servos  of 
future  gun  systems  will  be  developed  and  represented  by  these  tables  for 
evaluations  with  the  ISO-PK  program. 

Combining  Gun  System  Errors 

The  errors  from  the  gun  servo  system  are  assumed  to  add  to  the  gun 
pointing  errors  resulting  from  sensor  errors.  Therefore,  since  these 
two  error  sources  are  assumed  to  be  Independent,  their  statistics  can 
be  added  to  obtain  the  mean  and  covariance  of  the  total  gun  pointing 
error . 

The  gun  pointing  errors  produced  by  the  above  analysis  are  In  angular 
units.  However,  the  errors  required  by  the  BPK  subroutine  are  In  units 
of  length.  In  a coordinate  system  centered  at  the  center-of-vulnerablllty 
of  the  target.  This  coordinate  system  Is  perpendicular  to  the  line 
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from  the  gun  to  the  target  at  time  of  predicted  intercept.  The  conversion 
in  the  elevation  direction,  from  angular  units  to  meters,  is  achieved  by 
multiplying  the  elevation  error  by  the  range  (R)  to  the  target  at  pre- 
dieted  intercept.  In  the  azimuth  direction,  the  conversion  is  obtained 
by  multiplying  by  RCOS(E),  where  E is  the  elevation  of  the  target  at 
predicted  intercept.  The  cosine  term  is  needed  to  convert  the  azimuth 
error  in  the  ground  plane  to  an  azimuth-like  (transverse)  error  in  the 
elevated  coordinate  system. 

The  statistics  describing  the  errors,  in  meters,  in  the  elevated 
coordinate  system  are  as  follows: 

Among-Burst  Error: 


wiJVClL.  XdLi^C  • 

“ " 

a 2 COV(A,E) 

R COS(E)  0 

R COS(E)  0 

COV(A,E) 

0 R 

(COV(Z)  + COV(S)) 

0 R 

- 

• 

These  combined  statistics  are  those  required  to  describe  the  among-burst 
errors.  The  BPK  subroutine  uses  the  correlation  coefficient  instead  of 
the  covariance  term  COV(A,E).  This  correlation  coefficient  is  defined 
as: 


A,E 


COV(A,E) 

Ve 


• - 

Mean: 

- • 

R COS(E) 

R 

f 
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Wlthln-Burst  Errors 


A brief  derivation  of  the  within-burst  errors  is  presented  here, 
similar,  more  complete  derivation  can  be  found  in  reference  3.  The 
angular  residual  dispersion  error  and  range  variation  error  of  the  pro 
jectile  must  be  projected  into  the  elevated  coordinate  system  centered 
at  the  target.  Define  the  errors  as  follows: 


A.  - azimuth-like  (transverse)  angular  dispersion 
A 

Ag  - elevation  angular  dispersion 

A„  - projectile  range  error 
R 

These  errors  are  assumed  to  be  independent,  normally  distributed  with 

zero  means  and  variances  of  o^,  o^,  respectively.  The  variance  of 

a e r 

the  projectile  range  error  is  due  to  the  muzzle  velocity  variation. 


a 


2 

t 


(JL)2 

m 


a 


2 

V 


Where 

V * muzzle  velocity  of  the  projectile 
m 


“ variance  of  the  muzzle  velocity 

Since  there  is  an  error  in  the  range  coordinate,  the  time  to  equal 
range  of  the  target  and  projectile  must  be  determined.  This  time  can  be 
approximated  by  the  following  if  >>  R. 


Al 


AR 


Vo 


> 


3 

Banash,  Robert  C.,  Notes  of  the  University  of  Michigan  Analytic  Gun 
Model,  Technical  Note  SY-TN3-70,  Systems  Analysis  Directorate,  HQ, 
US  Army  Weapons  Command,  Rock  Island,  IL,  April  1970. 
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vrtiere 

R * range  rate  of  the  target 

V = velocity  of  the  projectile  at  range  R 
o 

The  velocity  components  of  the  target  in  the  elevated  coordinate 

• • • • % 

system  are  (ARCOS(E),  R,  ER)  * Where  A and  E are  the  azimuth  and 
elevation  rates  of  change  of  the  target  at  the  time  of  predicted  intercept. 
The  velocity  components  of  the  projectile  are  (0,  Vg,  0).  Therefore, 
the  relative  velocity  components  between  the  target  and  projectile  are 
-(ARCOS(E),  “Vg,  ER) . The  errors » in  meters,  at  the  time  when  the  pro- 
jectile and  target  are  at  equal  range  are 

AA'  = R AA  - At  ARCOS(E) 

= R AA  + ^ ARCOS(E) 

and 

AD  • 

AE ' * R AE  + ^ ER 
S 

AA'  is  the  error  in  the  traverse  direction  and  AE'  is  the  error  in  the 
elevation  direction.  Since  these  quantities  are  sums  of  independent, 
normal  random  variables,  they  have  a bivariate  normal  distribution  with 
the  following  parameters. 

Within-Burst  Errors: 

Means : 

E [AA'J  = 0 
E [AE']  = 0 
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Variances: 


CA'  - R' 


'e’ 


? . I RC0S(E)A  \ 2 a? 

5i  / 

v^l 


R^„l* 


Correlation  coefficient; 


^A’,  E' 


^ C0V(AA*.AEM 
^A’‘^E* 


'A'  E 


- — COV((R  AA  + (r  AR)) 

Vs  Vs 


‘^A'‘^E’ 


1 AE  R^  COS(E)  2 

V? 


At  this  point  all  the  necessary  inputs  are  presented  for  the  computa- 
tion of  burst-probability  of  kill  (BPK)  as  presented  in  the  form  in 
reference  2,  starting  at  page  6.  The  following  relates  the  notation  used 
in  this  report  to  that  of  reference  2: 

For  Within-Burst  Errors: 


Ml  = E[AA'] 
M2  = E[AE'] 


2 


Loc«  Cit. 
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For  Ainong-Burst  Errors: 


Before  presenting  the  equation  used  for  BPK  from  references  2,  3 or 
4,  the  basic  assumption  limiting  the  validity  of  this  equation  is 
reiterated: 

a.  Target  vulnerability  area  is  small  compared  to  the  within-burst 
dispersion. 

b.  Rounds  can  be  grouped  within  bursts  such  that  the  burst  centers 
can  be  considered  independently  distributed  while  the  rounds  are  distri- 
buted independently  about  the  burst  centers. 

The  ' notation  used  in  the  following  equations  means  transpose  of  the 
indicated  matrix. 

Under  these  assumptions,  the  equation  for  BPK  is  given  as: 


In  Air  Defense  Gun  Evaluation,  Technical  Report  No.  (CADES  3-1),  Artillery 
and  Air  Defense  Weapons  Directorate,  Rock  Island. Weapons  Laboratory, 

Rock  Island,  IL,  September  1971. 


BPK  - E C)  (-1) 

j-1  ^ 


n 


exp{-  ^ v'  A[l-j(I+jA)  ^ A]  v} 


where  n is  the  number  of  rounds  in  the  burst  and 


2 

Loc.  Cit. 

3 

Loc.  Cit. 

Banash,  Robert  C.,  A Markov  Chain  Approach  to  Modeling  Tracking  Error 
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A - K'P' 


P 

0}  1 2 


PK 


P is  the  matrix  of  characteristic  vectors  of 


V = 


P’VP 


Finally, 


V = k"^P 


0 


0 


M, 


M. 


= r2 


Ay  Is  the  target  vulnerable  area. 
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INPUT  DATA  DEFINITION 


The  input  data  is  divided  into  three  sections: 

(1)  data  to  control  the  planes,  (2)  gun  system  data,  at  which  the  plots 
are  to  be  made  and  (3)  target  data.  The  program  is  set  up  such  that 
the  data  for  each  category  is  read  from  a different  numbered  data  set. 

Each  category. is  explained  with  regard  to  data  content,  meaning,  and 
format , as  follows : 

Control  of  Selected  Planes. 

This  data  controls  the  size  and  location  of  the  planes  on  which  the 
footprint  is  printed.  It  also  determines  the  velocity  (magnitude  and 
direction)  of  the  target.  Each  card  is  described  in  the  order  in  which 
it  is  read.  In  the  following  discussion,  the  paragraphs  are  headed  by 
the  list  of  variables  whose  values  are  read.  Detailed  description  of 
the  variables  and  how  they  relate  to  the  program  are  also  presented.  This 
format  is  also  used  in  the  next  two  sections  — Gun  Systems  Data  and 
Target  Data. 

NVEL,  NCR,  NALT,  PLOTXY,  PLOTXZ 

(315, 2L5) 

NVEL  - The  number  of  different  velocities  which  will  appear  on  a 
later  card.  One  plot  will  be  made  for  each  velocity. 

NCR  - The  ntimber  of  different  crossing  ranges  at  which  XZ  plots 
are  to  be  made.  The  actual  values  of  the  crossing  ranges 
appears  on  a later  card.  Also,  this  data  only  has  meaning 
when  an  XZ  plot  is  to  be  made. 

NALT  - The  number  of  different  altitudes  at  which  X-Y  plots  are 
to  be  made.  Again  the  actual  values  of  the  altitudes 
appear  on  a later  data  card,  and  this  card  is  of  concern 
only  when  X-Y  plots  are  to  be  made. 
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PLOTXY  - This  logical  constant  should  be  set  equal  to  T (true)  if 
X-Y  plots  (at  various  altitudes)  are  to  be  made.  If  no 
X-Y  plots  are  required,  it  should  be  set  equal  to  F (false). 

PlOTXZ  - This  logical  constant  should  be  set  equal  to  T (true)  if 
XZ  plots  (at  various  crossing  ranges)  are  to  be  made.  If 
no  X-Z  plots  are  required,  it  should  be  set  equal  to  F 
(false) . 

XMIN,  XMAX,  YMIN,  YMAX,  ZMIN,  ZMAX 
(6F10.0) 

XMIN  - The  combination  of  all  the  data  on  this  card  sets  the  size 
of  the  plots  which  are  to  be  made.  XMIN  contains  the 
minimum  value  of  the  X-axis. 

XMAX  - The  maximum  value  of  the  X-axis. 

YMIN  - If  an  X-Y  plot  is  to  be  made,  the  minimum  value  of  Y (l.e., 
the  minimum  crossing  range)  is  contained  in  the  variable 
YMIN. 

YMAX  - The  maximum  value  of  the  Y-axis.  The  combination  of  XMIN, 
XMAX,  YMIN,  and  YMAX  describes  the  size  of  the  X-Y  plane. 

ZMIN  - If  an  X-Z  plot  is  to  be  made,  the  minimum  value  of  the 
Z-axls  (minimum  altitude)  is  contained  in  this  variable. 

ZMAX  - The  maximum  value  of  the  Z-axls.  The  combination  of  XMIN, 
XMAX,  ZMIN,  and  ZMAX  describes  the  size  of  the  X-Z  plane. 

VEL  (I),  I = 1,  NVEL 
(7F10.0) 

The  VEL  array  contains  the  actual  values  of  the  different  velocities 
for  which  plots  are  to  be  made.  The  number  of  velocities  which  will 
be  read  is  equal  to  NVEL. 

GRANGE (I),  I » 1,  NCR 
(7F10.0) 

The  GRANGE  array  contains  the  actual  values  of  the  different 
crossing  ranges  for  which  X-Z  plots  are  to  be  made.  The  number 
of  values  which  will  be  read  equals  NCR. 


22 


ALT(I),  I - 1,  NALT 

The  ALT  array  contains  the  actual  values  of  the  different  altitudes 
at  which  X-Y  plots  are  to  be  made.  The  number  of  values  which  will 
be  read  is  equal  to  NALT% 

Altogether,  there  will  be  (NVEL*(NCR  + NALT))  plots  made  if  both  PLOTXY 
and  PLOTXZ  are  initialized  to  T (true).  If  only  PLOTXY  is  equal  to  true, 
then  there  will  be  (NVEL*NALT)  plots.  Where,  if  only  PLOTXZ  is  equal  to 
true^  there  will  be  (NVEL*NCR)  plots.  The  size  of  all  plots  is  deter- 
mined by  the  parameters  on  the  second  card. 

Gun  Systems  Data 

The  data  contained  in  this  data  set  describe  the  parameters  related 

to  the  gun  system.  This  includes  the  gun  constraints  as  well  as  the 

description  of  the  random  errors  of  the  system.  The  same  format  will  be 

used  in  this  discussion  as  was  used  in  the  previous  section  in  describing 

the  data,  card  by  card. 

TITLE 

(20A4) 

This  is  an  80-character  verbal  title  which  appears  on  all  the  plots. 
It  may  contain  such  identifying  information  as  the  name  of  the 
system  being  Investigated  or  some  particular  configuration  of  a 
system.  The  actual  content  of  the  information  on  this  card  is 
optional  but  must  be  limited  to  80  characters. 

CALIB 

(15) 

The  value  of  variable  CALIB  is  used  to  determine  from  which  numbered 
data  set  the  target  data  is  to  be  read.  In  doing  this,  10  is  added 
to  the  value  of  CALIB  to  determine  the  data  set  number. 
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SIGVO,  SIGXD,  SIGYD,  VO 
(4F10.0) 


0 


The  standard  deviation  of  the  muzzle  velocity  of  the  projectile 
(o  ).  The  units  are  meters/sec. 


SIGXD 

- The  standard  deviation  of  the  residual  dispersion  of  the 
projectiles  in  the  traverse  direction  (^^)  • This  variable 
should  be  in  units  of  radians. 

SIGYD 

- The  standard  deviation  of  the  residual  dispersion  of  the 
projectiles  in  the  elevation  direction  (^g) • This  variable 
should  be  in  units  of  radians. 

VO 

- The  mean  value  of  the  muzzle  velocity  of  the  projectile  in 
meters  per  second. 

EDMAX,  ADMAX,  TREAC,  ENDS,  RMAX,  ELMAX 


(6F10.0) 

EDMAX 

- The  maximum  elevation  rate  which  can  be  achieved  by  the 
gun  in  radians  per  second. 

ADMAX 

- The  maximum  azimuth  rate  which  can  be  achieved  by  the  gun 
in  radians  per  second. 

TREAC 

- The  reaction  time  for  the  gun  system  in  seconds.  This 

reaction  time  is  defined  as  the  time  from  initial  detection 
of  a target  to  the  time  when  the  system  can  begin  firing. 

RNDS 

- The  number  of  rounds  in  a burst  from  the  gun. 

RMAX 

- The  maximum  effective  range  of  the  projectile  in  meters. 
This  is  usually  determined  by  the  range  at  which  the 
velocity  of  the  projectile  is  at  such  a low  level  that  it 
can  no  longer  inflict  damage  to  the  target. 

ELMAX 

- The  maximum  elevation  angle  in  radians,  to  which  the  gun 
can  elevate. 

R,  TF 
(2F10.0) 

These  two  variables  are  used  to  determine  the  ballistic  drag 
coefficient  of  the  projectile.  R is  a range  in  meters  and  TF 
is  the  time  of  flight  in  seconds  to  that  range.  These  two  para- 
meters can  be  taken  from  anywhere  on  the  range  - time  of  flight 
curve  of  a projectile,  especially  if  the  projectile  obeys  the 
"3/2  law"  for  the  drag  force.  If  the  projectile  doesn't  follow 
this  law,  a good  approximation  is  usually  achieved  at  a range 
which  is  2/3  the  maximum  range. 
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i. 


DYNAM,  DYNEM,  MAXAZ,  MAXEL 
(2F10. 0,215) 


DYNAM  - The  mean  gun  servo  error  in  the  azimuth  direction.  This 

error  is  usually  considered  to  be  a bore  sight  error  since 
servo  mechanisms  are  usually  not  designed  with  a mean  error 
in  them.  The  units  are  radians. 

DYNEM  - The  mean  gun  servo  error  in  the  elevation  direction. 

This  is  usually  the  result  of  a bore  sight  error.  The 
units  are  radians. 

MAXAZ  - The  number  of  entries  in  the  azimuth  rate  - standard 
deviation  tables.  (MAXAZ  ^ 10) 

MAXEL  - The  number  of  entries  in  the  elevation  rate  - standard 
deviation  tables.  (MAXEL  ^ 10) 

ELDOT(I),  I - 1,  MAXEL 
(7F10.0) 

The  ELDOT  array  contains  the  elevation  rate  entries  in  the 
elevation  rate  vs.  standard  deviation  tables.  The  units  are 
radians  per  second  and  the  number  of  entries  should  equal  the 
value  of  MAXEL. 

DYEL(I),  I - 1,  MAXEL 
(7F10.0) 

The  DYEL  array  contains  the  standard  deviations  (in  radians)  in  the 
elevation  direction  of  the  gun  servo  system  error.  The  quantities 
of  this  array  must  correspond  in  order  to  the  elevation  rates  in  the 
ELDOT  array.  The  combination  of  these  two  arrays  (ELDOT  - DYEL) 
describes  a curve  relating  the  standard  deviation  of  the  gun  servo 
error  in  the  elevation  direction  to  the  gun  elevation  rate.  A 
linear  interpolation  between  the  quantities  in  these  arrays  is  used 
to  achieve  the  standard  deviation  of  the  elevation  error  at  any 
particular  value  of  elevation  rate. 

AZDOT(I),  I - MAXAZ 
(7F10.0) 

The  AZDOT  array  contains  the  gun  azimuth  rate  entries  in  the  azimuth 
rate  vs.  standard  deviation  tables.  The  units  are  radians  per 
second  and  the  number  of  entries  should  equal  the  value  of  MAXAZ. 
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DYAZ(I) , I - MAXAZ 
(7F10.0) 

The  DYAZ  contains  the  standard  deviations  (in  radians)  in  the 
azimuth  direction  of  the  gun  servo  system  error.  As  was  the  case 
with  the  DYEL  array,  the  entries  in  the  DYAZ  array  should  correspond 
to  the  rates  in  the  AZDOT  array,  and  a linear  interpolation  is  used 
to  find  the  standard  deviation  at  any  particular  value  of  azimuth 
rate. 


XM  (I) 
(6F10.0) 


The  XM  array  contains  the  mean  value  of  the  measurement  errors. 
Earlier  in  the  text  of  this  report  these  quantities  were  described 
with  the  symbol  IT  . The  units  are  radians  - meters  - seconds. 

The  order  is:  1)  range,  2)  azimuth,  3)  elevation,  4)  range 

rate,  5)  azimuth  rate,  6)  elevation  rate. 

SIG  (I,  J) 

(6F10.0) 


The  SIG  matrix  contains  the  standard  deviations  and  correlation 
coefficients  of  the  measurement  errors.  The  terms  on  the  diagonal 
of  this  6x6  symmetric  matrix  are  the  standard  deviations  of  the 
errors.  The  off  diagonal  terms  are  the  correlation  coefficients. 
The  matrix  is  read  in  one  row  at  a time  (l.e.,  one  row  per  card) 
following  the  same  order  as  the  mean  errors.  This  is,  range  first, 
azimuth  second,  etc.  The  units  are  radians  - meters  - seconds.  To 
Illustrate  the  order  of  the  cards,  the  information  on  the  first 
three  cards  is  as  follows: 


1st  card  - a 


r,a  r,e  r,r  r,a 


r,e 


2nd  card  - p cr  p 


P-  : P- 


r,a  a a,e  a,r  a, a a,e 


3rd  card  -p  p ap*p*p* 
r,e  a,e  e e,r  e,a  e,e 


etc. 


Target  Data 


The  target  is  described  as  a combination  parallelepipeds  which 
represent  the  vulnerable  components  which  make  up  the  target.  There  may 
be  as  many  as  seven  components  which  make  up  the  target.  Three  values 
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of  each  of  the  six  faces  of  the  components  are  input  to  the  program.  The 
three  values  of  each  area  correspond  to  the  vulnerable  area  presented  as 
a function  of  projectile  striking  velocity.  Presently,  the  program  is 
set  up  so  that  the  first  area  for  each  face  is  for  striking  velocity 
of  500  ft/sec  and  the  second  area  is  for  1000  ft/sec.  The  third  value  is  not 
currently  being  used;  however,  space  must  be  allowed  for  it  on  the  data 
cards.  As  was  mentioned  in  the  previous  section,  the  value  of  the 
variable  CALIB  is  used  to  determine  the  numbered  data  set  on  which  this 
target  data  is  found.  A card  by  card  description  of  the  data  is  given 
below. 

EMAX 

(15) 


The  number  of  vulnerable  components  from  which  the  target  is  made 
should  be  read  into  the  variable  EMAX. 


S(I,  J,  K) 

(6F10.5) 

The  S array  contains  the  values  of  the  areas  of  the  faces  of  the 
vulnerable  components  in  units  of  square  meters.  Because  of  the 
format  by  which  these  quantities  are  read,  three  cards  are  required 
to  define  each  of  the  (EMAX)  components.  The  first  three  quantities 
on  the  first  card  are  the  three  areas  of  the  front  side  of  the 
component , while  the  second  three  are  the  areas  of  the  rear  of  the 
component.  On  the  second  card  are  the  areas  of  the  port  and  star- 
board sides  of  the  component  (in  that  order) . And  on  the  last  card 
the  areas  of  the  bottom  and  top  sides  of  the  component  are  given 
(in  that  order). 


Next  page  is  blank. 
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LOGIC  DIAGRAM 


A schematic  diagram  of  the  program  and  subprogram  logic  is  presented 
in  Figure  5. 


Next  page  is  blank. 
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MAIN  PROGRAM 


Figure  5 


Logic  Diagram  (1  of  8) 


Figure  5.  Logic  Diagram  (2  of  8) 
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SUBROUTINE  PLOT 


ft 


Initialize  X-Axis 
And  ”B"-Axls  For 
Plotter  Subroutine. 
Initialize  Plotting 
Matrix. 

1 

Initialize  Save-Matrix 
For  Burst  Kill  Prob. 
(BPK) 


1 

Initialize  X to  Xmln 
For  Sweep  Left  To 
Right 

I 


( 


Do  Thru  a 
For  IX-1,  101 


T 


} 


Compute  BPK  At 
One  Step  Down  From 
Bmin.  Set  to  PKS 


Initialize  "B”  Variable 
To  Bmin  For  Sweep  From 
Bottom  To  Top 


Do  Thru  6 
For  IB-1,  51 


i 

Call  COMPK  To 
Compute  BPK  At 
Current  Grid  Point 
And  Save  As  PKC 

i 

”SaveTK^TfiT""”" 

Matrix  PSAVE 


5 


Figure  5.  Logic  Diagram  (3  of  8) 
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^ Is  Slant  Range  At  Pre- 
I dieted  Intercept  Greater  k 

^han  Max.  Gun  Range?  T 

I No 


Yes 


1 


Compute  Range  At  Predicted 
Intercept  After  The  Target 
Has  Flovm  ”TREAC"  Seconds 
At  Same  Vel. 


c 


I 


(First  GrT^Tin^  „ 
On  Plot? 


Plot  ”R*' 


Is  New  Range  Greater 
Than  Max.  Gun  Range? 


> 


Yes 


T 


1 


No 

Plot  "T”  In 
Current  Grid 
Position 

. 1 

Yes 

^ Is  Azimuth  Rate  GreateX 

1.  ildii  XXclAtt 

No 

X 

f 

Plot  "A’*  1 

Is  Elevation  Rate  Or 
Elevation  Greater  Than 
Max.  For  Gun? 


) 


Yes 


( 

c 


T 


1 No 

Set  SS  Equal  To  Sign 
Of  BPK  Gradient 

Plot  **E" 

5^ 


Do  Thru  Y 
For  IP=1,  NPK 


I 


Did  PKS  Just  Cross 
The  PK(IP)  Contour 
Value? 


) 

> 


Yes 


Plot  Corresponding 
Contour  Letter 


(f) 


Figure  5.  Logic  Diagram  (4  of  8) 
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Figure  5.  Logic  Diagram  (5  of  8) 
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SUBROUTINE  COMPKl 


c 


Read  Description  Card 
And  All  Gun  System  Data 






Write  Out  All  Input  Data 
With  Format  Titles 

Call  Subroutine  DYNGNl 
To  Read  In  Gun  Servo 
Errors. 

Call  Subroutine  MESIGl 
To  Read  In  Measurement 
Errors. 

Call  Subroutine  VULAPD 
To  Initialize  Vulnerability 
Data  In  VULNAR 

] 

i 

Compute  Squ; 
Used  In  Coi 

ared  Values 
input  ing  BPK 

Compute  Ground 
Range, 

Range  Squared, 
Range,  Elevation, 
Azimuth,  Range  Rate, 
Azimuth  Rate  And 
Elevation  Rate 


Call  Subroutine  GFUNC 
To  Calculate  The  Target 
Position  At  Predicted  Inter- 
cept And  the  Necessary 
Partial  Derivatives  Of 
The  Gun  Position 
Functions. 


c 


Is  Range  At  Predicted 
Intercept  Greater  Than 
Max.  Gun  Range? 


JL  No 

o 


> 


yes 


Set  BPK=0 


Return 


Figure  5.  Logic  Diagram  (6  of  8) 


Perform  Table  Look-Up 
For  Dynamic  Gun-Pointing 
Error  Standard  Deviations 
In  Azimuth  And  Elevation 


Compute  Bi 
Azimuth-L 
Elevation-L 
For  B 

ases  For 
ike  And 
ike  Plane 
URPK 

j 

r 

Compute  C 
Standard  D 
For  B 

ombined 

eviations 

URPK 

Call  Subroutine  UULNAR 
For  Total  Projected 
Area  Of  Target 


Call  Subroutine  BURPK 
For  Burst  Kill 
Probability  BPK 


Figure  5 


Logic  Diagram  (7  of  8) 
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Figure  5.  Logic  Diagram  (8  of  8) 
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FORTRAN  PROGRAM  SOURCE  LISTING 


This  section  provides  a computer  printout  of  the  FORTRAN  statements 
making  up  the  ISO-PK  program. 


Next  page  Is  blank. 
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o‘>or»oooooooooooooooooorioooooor>oooonoooooooooooooooooo 


FOPTRAN  IV  G LEVEL  21 


MAIN 


DATE  * 74109 


IS/IS/S*^ 


unimap  --AIR  DEFENSE  FOOTPRINTS 
THIS  PROGRAM  DEVELOPS  LINES  OF  CONSTANT  RURST-KILL  PROHABILITY  M)R 
AIR  DEFENSE  GUN  SYSTFMS.  THE  MODEL  USED  IS  AN  EXPANDED  VERSTOn 
THE  UNIVERSITY  OF  MICHIGAN  ANALYTIC  MODEL  OEVELOPFO  iN  THE  aFaaOS 
GUN  STUDY 

THE  main  program  acts  AS  A DRIVER  AND  CONTROL  FUR  THF  MODULES  Ihm 
PLOT  ANO  COMPUTE. 

THE  PROGRAM  READS  PLOTTING  CONTROL  INFORMATION  aNO  OPTIONS 
AND  PRINTS  TABLES  OF  EDUIVALENTS  BETWEEN  PROP AH  I L I T TFS  AND 
PLOTTING  CHaRACTFRS. 


coordinate  system  FOR  UNIMAP 
♦ 71 

i 

♦ Y 

/ 

/ 


-X 


/ #<-...  Target  posittdn 

/ / 

/ / 

/ / 

/ / 

/ / 

/ / 

/ / 

/ / / 

/ / / 

/ / / 


GUN 


azimuth  angle  measured  from  *X  IN  XY  PLA»^'E. 

elevation  angle  MFASURFO  from  XY  plamF  Tn 
SLANT  range  line  FROM  GUN  TO  TARGET. 


-Y/ 

/ 


#«###«•##««##««*««#««#« 


✓ 


consider  the  XY  PLANE: 


♦ Y 

) 

I 

I 

I 

I < < < 

I 


A1 


FORTRAN  TV  G LEVEL  21 


MAIN 


DATE  = 74109 


OftOl 

000? 

0o03 


0004 

OoO«S 

OOOG 

0007 

OOOR 

OOOQ 

OfllO 

0011 

0012 

oon 

0014 


c 

1 

c 

1 

< — 

target 

c 

1 

PA  THS 

c 

1 

r. 

1 

c 

1 

c 

1 

c 

1 

c 

1 

c 

1 

c 

1 

c 

1 

c 

C -X  (0»0>  ♦» 

C GUN 

c The  scenario  in  the  x?  plane  is  similar  but  with  i hf. ing  the 
c vertical  axis. 
c 

c 

common  X.  Y ( 1 I .zn  ) .XDOT.YOOT.ZDnTtXDDtynDtZDOtt  GDOI , AGOOT  tW'-.Uf  Vi>. 
. TF.RNOS. VO.RMSQf TREACtEOMAX. AOMAXtXMIN.XMAX, T ITLE  <20>  tCALl'* 
COMMON  /SPLOT/  PKtCHARf MCHAR*NPk 
INTEGER  CALIB 


• ■C 


c 

CALIH=1 ! 

?OMM 

HF 

PHOJ. . 

KK  KILL. 

c 

CALIH=?: 

30MM 

HE 

PROJ.f 

KK  KILL. 

c 

C4LI8=3: 

4 0 MM 

HF 

PROJ. • 

KK  KILL. 

c 

CALIH=4: 

S7MM 

HE 

PROJ. f 

KK  KILL. 

c 

CALIBsS: 

n LP 

WAPHEAOf 

KK  kill. 

c 

calib=g: 

?0MM 

HF 

PROJ. f 

CLEAN  aircraft,  a kill. 

c 

CALIHsT: 

20  MM 

He 

PROJ. f 

aircraft  with  wing  tanks. 

A 

c 

CALIBrfl: 

lOMM 

HF 

PROJ, • 

CLEAN  aircraft.  A KILL. 

c 

CALIH=9» 

30MM 

HF 

PROJ. f 

aircraft  with  wing  tanks. 

A 

c 

CALlHslO; 

40MM 

HF 

PROJ. 

. CLEAN  aircraft,  a kill. 

c 

CALIH=11 : 

40MM 

HF 

PHOJ. 

. aircraft  with  wing  Tanks 

c 

CALIB=12: 

5/MM 

HE 

PROJ. 

. CLEAN  aircraft.  A KILL. 

c 

CALIB=13; 

57MM 

HF 

PROJ. 

, AIRCRAFT  with  wing  TANKS 

f 1 

c 

CALI8=14: 

total  projected  AREA  OF  A FOREIGN  AIRCRAF 

c 

CALIBslS: 

GLaaOS 

TARGET 

: 4 so  M FRONT  & HACK!  ?4 

SQ 

KILL. 

KILL. 

KILL. 

KILL. 

(FOR  HIT  PROM 
SIOES.  TOR  X 


11 


real  VEL(IO) , alt (10) .CHANGE (10) .PK(IO) 
integer  *2  CHAH(IO) .MCHaR(4) 

LOGICAL  RLOTXY.PLOTXZ 

REAO  (5. 1 1 .EN0=1?)  NVEL.NCR.NALT. RLOTXY.PLOTXZ. 

XMIN.XMAX. YMIN. YMAX.ZMIN.ZMAX 
FORMAT  (3Ib,2Lb/(6F10.0) ) 

READ  (5.2»END=1?)  ( VfL ( I ) . 1=1 .NVEL) 

(CRANGE(I) .1=1 .NCR) 

(ALT(I> ,I=1.NALT) 


REAP  (5.?.EN0=12) 
REAO  <5.?.END=1?) 
? format  (7F10.0) 

13  WRITE  (1  »H 
1 FORMAT  CIG 


A n F S ISO 


PK  CONTOUR  PROGRa 


C'LC. 
• •or  Him 


42 


1 


forthan  tv 


OOl*; 

OOI<S 

0017 

ooie 

0019 


oo?o 

0021 

0022 

0023 

0024 
002«; 
002ft 

0027 

0020 

0029 

0030 

0031 

0032 


0033 

0034 

0035 
003ft 
0037 
003R 

0039 

0040 

0041 

0042 


0 LFVliL  21  main  DA7F  =s  74109  Ib/lft/SS 

march  1973') 

C CALL  SURROUTINE  TO  READ  PARTICULAR  GUN  SYSTEM  DaTA. 

CALL  COMPKl  (U2) 

WRITE  (l*fl>  (I»PK(n tCHAR(l) «I=1*NPK) 
ft  format  (•l*//»7X.'C0NT0UR*tftX.*PK*t7X.*CHAPACTLR*/ 

.('  '.IlO.FlB.A.RXAin 

WRITF  (1*10)  MCHAH 

10  FORMAT  ( • 0* « 1X« ‘MASK* . 19X, 'CHARACTER*/'  EL.  RATE  OR  max  FL'.T30.A1 
./•  AZ.  rate* «T30. Al/'  INSUF,  T I ME • ♦ T30 1 A 1 / ' OUT  OF  R ANGK • » T TO . A 1) 

C start  no  LOOPS  for  plotting  with  parametfrs  of  altitude. 
c target  velocity  and  crossing  range. 
no  3 IVEL=1'NVFL 
XOOTsVELdVEL) 

YOOT=0,0 
Z0OT=0.0 
XODzO . 

YDD»0. 

Z00*0. 

C TEST  IF  MO  "GROIJNO  PLOT"  OFSIREn. 

IF  ( .NOT.PLOTXY)  GO  TO  24 
C PLOT  FOR  "NALT"  target  ALTITUDES. 

DO  2S  IALT=1,NALT 
Z(1)=ALT(IALT) 

CALL  PLOT ( YMIN.YmaX.Y. IftHCROSSING  RNG.t  M) 

2S  WRITE  (1.7)  title. XOOT.Zd) 

7 FORMAT  ( *0' .20A4/ 

.*  target  VEL0CITY=* .FH.3. *M/S* .SX'TARGET  ALTI TUDE=* .F9.3. *M' ) 

C lEST  IF  *10  "SKY  PLOT"  DESIRED. 

24  IF  ( .NOT.PLOTXZ)  GO  TO  3 

C plot  in  vertical  plane  with  "NCR"  target  crossing  ranges. 

DO  27  ICP=1.NCR 
Yd)  =CRAMGE  (ICR) 

CALL  PLOT (ZMIN.ZMaX.2, IftHALTITUOE.  M£TERS) 

27  WRITF  (1.26)  title. XOOT.Y(l) 

2ft  format  (*0*.20A4/ 

.'  Target  velocity=*,f8.3. 'm/s* .sx'crossing  Range=* .f9.3. *m* > 

3 CONTINUE 
GO  TO  13 

12  call  exit 

END 
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fortran  TV 


0001 


Ono? 

0003 

0004 

0005 

0006 


0 LFVEL  21  main  rate  = 7410R  15/15/55 

C 

HLOCK  DATA 
C 

C I''J1TIALIZE  PLOTTING  CHARACTERS  AND  RK  CONTOUR  VALUE';. 

COMMON  /$PLOT/  PK fCHARtMCHAW.NPK 
INTFRFR  NPK/10/ 

REAL  PK( 10) /n. 005. 0.01 tO. 05*0.1 0«0.?5*0.5.0,7bf 0.9» 0,95. n.9<}/ 
INTFGFR  *2  CHARI  10) / 'F' . 'G* . »H' ♦ • I • ♦ • J* » «K» * 'L' » 'M» , IN' . 'O'/. 

, MCHaR<4)/*F'»*A».»T'.*R»/ 

END 
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ropTRAN 


0001 


Ooo? 

onn.i 

0004 

Opor 

Onn^ 

0O07 

OoOfl 

OOOP 

0010 

0011 

001? 

0013 

0014 

0015 
OOl  ft 

0017 

0018 
OOlP 

00?n 

00?1 

00?? 

00?3 

00?4 

00?5 

00?ft 

00?7 

00?P 

00?P 

0030 

0031 
003? 

0033 

0034 
0038 
0P3ft 

0037 

0038 

0039 

0040 

0041 


TV  6 l.EVEL  21  main  DATE  * T4109  15/15/8S 

r 

SUPROUriNE  PIJRPK  (AM3«XM4tSl  «S?  .S3  »S4  . R . XN . A V . RKPS ) 

C 

C 

C THIS  SUHROUTINE  COMPUTES  HURST  KILL  PROHARILITT.  ftPK,  USING  THE 

C analytic  model  started  HY  THE  UNIVERSITY  OF  MICHIGAN  AND 

c completeo  hy  this  command. 

C reference:  "NOTES  ON  THF  UNIVERSITY  OF  MICHIGAN  ANALYTIC  GUN 

C MODEL".  TECH  NOTE  SY-TN3-70.  ROBFRT  C.  BANASH.  SYSTEMS  ANALYSIS 

C DIRECTORATE.  HO.  USA  wECOm,  aPRR.  1970. 

C 

IMPLICIT  REAL*fl( A-H.O-Z) 

real  XM3.XM4.Sl .S?.S3.S4  »R. XN. AV.BKPS.X (?) 
peal  S0HI.AMIN1.AHS 

REAL*H  IIAl.JlA?.IIA3,IIAAl.IIAA2.irAA3 
DIMENSION  ALF(?)  .P(?.2)  .BKPJOOO) 

IF(  AHS(R) .LE.0.05)  GO  TO  2 
DO  1 J=l.2 

ALF(J)  =.S*(  S3+S4* (3-?*J) * SORT( (S3*S4)**2-4.*(i.-R**?)*S3*S4) ) 

X(J)  xR*  SORT  (S3*S4)/(aLP( J)-S3) 

P(l.J)  X S0RT<l./(l./X(J)**2»l.)) 

P(?.J)  = P(1.J)/X(J> 

I aLE(J)  X 1/aLF(J) 

GO  TO  3 

? ALF(l)  X 1./S3 
ALF(2)  = 1./S4 
P(l.  1 ) = 1. 

P< 1 .2)  X 0. 

R(2. 1 ) X 0. 

P(?.2)  = 1. 

3 RKP  X 0.0 
N X XW 

K X 1 
EX  X 0. 

C = 1.0 
FJMOal .0 
FNPOxN.l 

PJ  = AV/(ft.2831fl5*SORI(Sl*S2)  ) 

IF(PJ.GT. 1.000)  GO  TO  100 

IF  ( XM3.EQ. 0. . and. XM4.E0. 0. ) K x 2 

AlNx  (P(1,1)**2/S1  ♦ P(2,l)**2/S2  ) / ALF(1) 

A2Nx(P(l,l)*P(l,2)/Sl*P(2.2)*P(2.1)/S2)/0S0Rr(ALF(l)*ALF(2) ) 

A3Nx  (P(l .2)»*2/Sl  ♦ P(?.2)**2/S?  )/ALF(2) 

A2N?=A2N*A2N 
no  29  Jxl.N 

NNxj 
FJ  X J 
FJlxl.n/FJ 

nENOTAxDSORT( ( A 1 N*F J1 ) * ( A 3N. F J1 ) -a2N2) 

GO  TO  (4,5) .K 

4 Al  X A1N*FJ 
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1 


FOpTRAN  tv  G LFVEL 


OOA? 

00A3 

0044 

004G 

004G 

0047 

004P 

004Q 

OOGO 

0051 

OOS? 

OOSl 

00G4 

00G5 

OO^^G 

0057 

005R 

00*^9 

OOfO 

00^1 

00ft? 

0053 

0054 

0055 

0056 

0057 
005B 
005Q 

0070 

0071 
007? 


21 

A2  = 

A2N*FJ 

A3  = 

A3N*FJ 

DIA  ^ 

= (Al*l 

IlAl 

= (A3* 

miWPK 


DATE  = 7410^ 


?9 

11 


1? 


100 


/ DIA 

XIA?  = -A?/DIA 
IIA3  = (Al*l*)/flIA 
IIAAl  = lIAl^Al^T IA2*A? 

IIAA?  = IIAl*A2^IT42*A3 
11AA3  = lIA2*A2«IIA3*a3 

API  = Al* ( I .“I  I aAl ) ~A2*T T aA? 

AP2  s -a1 *11 AA2*A2* ( 1 .-T I aa3) 

AP3  = -A2*IIAA2*a3»(l.-ITAA3) 

Vl=0S0RT(ALP(in*(P(l*l»*XM3*P  2.1  *XM4 

\/?-n<;ORr  (ALF  (2)  ) * (P(  1 f2)  *XM3  + P(2»2) 

FX~=-(VI*V1*AP1*2*V2*V1*AP2*V2*V?*AP3)*0.‘> 

TF(1)A«S(FX)  .GT.170.000)  EX=0SIGN(  170.000. EX) 
C*-(  (FNP0-FJ)/FJ)*PJ*(F.IP0/FJ)  *c 
BKP,I(  J)  =r*OEXP  (EX)  / DENOIA 

XF  (naBG(BKPJ(J) ) .LE. 1 .00-10)  GO  TO  11 

CONTTNUF 

CONTINUE 

no  12  J=1.NN 

PKP=HKP-PKPJ(NN-J+ 1 ) 

CONTINUE 

HKPS  = DMlNl (l.ODO.BKP) 

PETUPN 
CONTINUE 
RKPS=1 .0 
HETUPN 

ENn 


15/m/GG 
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FORTPAM  TV  R LFVFL  21 


MAIN 


OATF 


74109 


IS/IS/S*^ 


0001 


FUNCTION  TLUl (X.YtN.XO) 


THI<;  SUHPOUTINF  PERFORM*;  TAHLE  LOOK-UP  AND  LINEAR  If4TFRPOLAT  ION 

ON  sinoly-surscpipted  arrays. 


000? 

oimfnsion  x(n) .y(n) 

0003 

IF  (N.LF.l)  00  TO  1 

0004 

IF  ( XO.LF .X { 1 ) ) 00  TO 

ooos 

DO  ? I=?.N 

OoOAv 

IF  (XO.LF. X(D)  00  TO 

0007 

? 

CONTINUE 

OOOfl 

I = N 

OOOP 

.1 

TLUl=Y(I-n*(XO-X(I-l 

0010 

return 

0011 

1 

TLUlsY(l) 

001? 

RETURN 

0013 

END 
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FORTRAN  IV  G LEVEL  21 


MAIN 


DATE  * 74109 


ooni 


0n02 

0n03 

0004 


0005 

OnOG 

Ono7 

Oooo 

0009 

00)0 

0011 


0012 

0013 


0014 

0015 


OOlG 

0017 


Oolfl 

O019 

00?0 

0021 

0022 

0023 


SUPROUTINE  OYNGNl 

r 

common  niJMY  (9)  »EGDOT.4GnnT 

COMMON/D YNERR/  OYNAM , OYNEM t DYNAV » DYMFV 

niMFNSION  OYAZ(IO)  tAZOOTOO)  tOYEL(lO)  .ELDOT(IO) 

C READ  IN  DYNAMIC  GIJN  POINTING  POSN.  ERROR  INFORMATION 
C read  in  tables  FOR  GUN  POINTING  ERRORS  IN  AZ.  AND  EL. 

C 

RE AD (2. 3)  DYNAM.OYNEMf MAXAZ.MAXFL 
3 FORmaT(2F10. 0.215) 

read  (2,2)  (ELOOT ( I) . 1=1 .MAXED 
read  (2,2)  (DYFL(I) .1=1, MAXED 
REAO  (2,2)  ( AZnoT ( I) , 1=1 .MAXAZ) 

REAO  (2.2)  (DYA2 ( I ) . 1=1 .MAXAZ) 

2 format  (7F10.0) 

C 

C WRITE  OUT  ALL  INPUT  DATA 

C 

WRITE  (1.20)  DYMFM,  (l.ELDOTd)  .OYEL(I)  .1  = 1. MAIL'D 

20  format  CODYNAMIC  GUN-POINTING  ERROR  FUNCTION  - ELEVATION* 
.5X*BIAS.  HAD.=*.F7.4/ 

.*  POINT* .2X13HEL.  RATE.  P/S5X23HSTO  OEV  OF  EL.  ERROR.  »/ 
.(*  *.I5,F1!>.3.F20.4)  ) 

WRITE  (1.21)  OYNAM, (I.AZOOT ( I ) ,OYAZ (I ) . 1=1 .MAXAZ) 

21  format  (*OOYNAM1C  GUN-POINTING  ERROR  FUNCTION  - AZIMUTH* , 
.7X*BIAS.  RA0,='.F7.4/ 

.*  POINT  * .2X13HA7.  RATE.  R/S5X23HSTD  DEV  OF  AZ.  ERROR,  H/ 
.(*  * .I5.ri5.3,F20.4) ) 

RETURN 

ENTRY  OYNGUN 

C CALCULATE  MEAN  AflO  VARIANCE  OF  DYNAMIC  POINTING  POSN. 

C 

C PERFORM  table-look-up  TO  OBTAIN  STANDARD  DEVIATIONS  OF  GUN 

C POINTING  ERRORS  IN  AZIMUTH  AND  ELEVATION. 

C 

5IGDYA=TLU1 ( AZDOT , DY AZ . M AX AZ , ABS ( AGOOT ) ) 

OYNaV  = SIGDYA**2 

SIG0YE=TLU1 (FLDOT,DYEL.MAXEL.ABS(E6DOT) ) 

0YNFV  = SIGDYE***2 
RETURN 

end 


15/15/5K 


48 


FORTRAN  tv  R LEVFL  21 


MAIN 


RATE 


^4109 


1 b/lS/=.S 


r 

n 

0001 

L 

SUHROIITINE  TFLlfTF.*) 

r 

rxirrr  3:=  = = = = x ===r  = = ===  = = « = r = = = = = = = = = = ===  = = = = = T = = = = ==  ==-  = 

000? 

c 

C0MM0N/Z7TFLT/  C0F(4) 

0003 

c 

dimension  R(4) 

c 

THIS  ROUTINE  SOLVES  THE  ANALYTICAL  TIME  OF  FLIGHT  EOllATION 

c 

THE  COEFFICIENTS  i)F  THE  FOURTH  OROER  TIME  OF 

c 

flight  FOUATION  ARE  STORED  IN  COF(I) 

r. 

SOLVF  the  FOURTH  ORDER  FOR  ITS  HOOTS 

0004 

CALL  URTIClCOFfRfNRE) 

L 

c 

FIND  THE  smallest  REAL  ROOT 

onn*; 

L 

TF  B 1.0F06 

L 

c 

ARE  there  FOUR  REAL  ROOTS  ?7? 

OOOR 

r 

c 

IFINRE  ,FQ.  4)  60  TO  5 

ARE  there  no  real  hoots  ??? 

c 

Ono7 

c 

c 

lElNRE  .FO.  0)  GO  TO  14 

THEHE  are  only  two  heal  ROOTS 

0008 

L 

GO  TO  o 

0009 

s 

IF(H(4))  7»7t6 

0010 

6 

TF  = AMim  (TF«H(4)  ) 

0011 

7 

iF(Hn))  9t9.e 

001? 

a 

TF  a AMINl (TF,R(3) ) 

O013 

<? 

IF(R(?n  lltlltlO 

0014 

10 

TF  a AMINl (TF.H(2) ) 

0019 

1 1 

IF(Rtn)  13»13.1? 

0016 

1? 

TF  a AMINl (TFtR( 1 ) ) 

0017 

13 

IF(TF  .OF.  1.0E06)  GO  TO  14 

0018 

RETURN 

L 

C 

p 

ThEH  are  NO  REAL  ROOTS  OR  THERE  ARE  NO  POSITIVE  REAL  ROOTS 

0019 

U 

RETURN  1 

oo?o 

END 
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FOPTRAN  IV  fi  level  21 


MAIN 


RATE  = 74109 


15/lS/S 


0001 


0002 

0003 

0004 

0005 

0006 
0007 
OoOfi 

0009 

0010 
0011 
001? 

0013 

0014 

0015 

0016 
0017 
0ol« 

0019 

00?0 

0021 

00?2 

00?3 

00?4 

OOPS 

0026 

00?7 

002fl 

0029 

0030 

0031 

0032 

0033 

0034 

0035 

0036 

0037 
0036 
0039 


C 

SUMROUTINE  QRTICtCtRfNRE) 

C 

C SOLVES  A polynomial  eouation  of  the  type 
c x**4  ♦ cm*x**3  ♦ C(?)*x**2  ♦ C(3)*X  ♦ C(4)  = 0 

C THE  COEFFICIENT  OF  X**4  IS  ASSUMED  TO  »E  1 

C R CONTAINS  THE  ROOTS 

C NRE  CONTAINS  IHE  NUMBER  OF  REAL  ROOTS 

C IF  THERE  ARE  TWO  REAL  ROOTS  THEY  WILL  BE  IN  R(l)  AND  R(2) 

C WITH  THE  COMPLEX  ROOTS  R(3)  ♦-  R(4)*I 

C IF  THEPE  ARE  NO  REAL  ROOTS*  THE  COMPLEX  ROOTS  ARE 

C R(l)  ♦-  R(2)*I  AMD  R(3)  ♦-  R(4)*I 

C 

DIMENSION  C(4) *R<4) *CP(3) .Y(3) 

ClS(J=Ctl)**2 

CPm=-C(2) 

CP(?)=C(1)*C(3)-4,*C(4) 

CP(3)=14.*C(2)-CIS0)*C(4)-C(3)**2 

call  cubic(cp*y,nre) 

A=ClSO/4.-C(2)*Y(l) 

B=.5*C( 1)*Y(1)-C(3> 

D=.?5*Y(1)**2-C(4) 

IFIA.RT.O.)  60  TO  10 
F=0. 

GO  TO  20 
10  E=SQRT(A) 

20  IF(D.GT,n.)  GO  TO  30 
F = 0. 

GO  TO  50 

30  F=SI6N(S0RT (0) *H) 

50  NRE=0 

REAI.a-,25*C(l)*.S*E 
0SCR=REAL**2-.5*Y ( 1 ) «F 
RAOsSORT (ABS(OSCR) ) 

IFIOSCR.LT.O. ) GO  TO  60 
NRE=? 

R(1)=REAL*RAD 
R(2)=PEAL“HA0 
GO  TO  6S 
60  R(3)=PFAL 
R(4)=RAD 
65  REAL=REAL-E 

D5>CR=REAL**2-.5*Y  ( 1 ) -F 
RAD=SORT(ABS(OSCW) ) 

IFCOSCR.LT.O.)  go  to  bo 

NREaNRE+2 

R(NPF»=E 

R(NRE)=REAL-RAD 

R(NRE-1)=REAL»RA0 

RETURN 

80  R(NPE*1)=REAL 
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R(NRF*?)xRAD 

RETURN 

end 
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DATE 
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0001 


0002 

Ooni 

0004 

0005 
OoO^< 
0007 
OoOft 

0009 

0010 
0011 
0012 

0013 

0014 

0015 

0016 
0017 

ooia 

0019 

0020 
0021 
0022 

0023 

0024 

0025 

0026 
0027 
002R 


C 

subroutine  CUflIC(C»R»NRn 

c 

C SOLVES  A polynomial  EQUATION  OF  THE  TYPE 

C X»*3  ♦ C(1)*X**2  ♦ C(2)*X  ♦ C (3)  = 0 

C THE  coefficient  OF  A**3  IS  ASSUMED  TO  BE  1 

C P CONTAINS  THE  ROOTS 

C NRE  contains  the  NUMBER  OF  REAL  ROOTS 

C IF  THERE  IS  ONE  REAL  ROOT  IT  WILL  RE  R(l)» 

C with  the  COMOlex  roots  R(2)  ♦-R(3)*I 

C 

niMENSION  C(3) .R(3) 

C1S0=C(1)**2 

P*C(2)-ClSU/3. 

Q=C(3)-C(1)*(C(2)/3.-2,*C1S'1/27.) 

DEL=4.«P*‘»3*27,*0**2 

T*C(1 )/3. 

lECOEL.LT.O.)  GO  TO  10 
SO=SORT<DEL/10B.) 

HQ=.5«0 

As-HO*SO 

B=-HQ-S0 

CRTa=SIGN(AHS(A)**(1.0/3.0) .A) 
CRTR=SIGN(ABS(R)**(1. 0/3.0) fH) 

Y*CRTA+CRTB 

R(1)=Y-T 

R(2)=-.5*Y-T 

R (3) =,066025404* (CRTA-CPTR) 

NRE=1 

RETURN 

10  PHI3=ATAN2(S0RT(-DEL/27.) .-0)/3. 
CON=?.*SORT (-P/3. ) 

R(l)=CON*COS(PHI3)-r 
R (2) =-COM*COS ( 1 . 04719B-PHT3) -T 
R(3)s-C0N*C0S(1 ,04719R*PHI3)-T 
NRE  = 3 

return 

END 
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RATE  = 74109 


0001 


OnO? 


0003 

0004 
0009 
OOOf. 
0007 
OoOfl 

0009 

0010 


0011 

001? 

0013 

0014 


0019 

OOliS 

0017 

OOlfl 

0019 

0020 

0021 

002? 

0023 


C 

C 

C 

C 

C 

C 

C 


C 

C 

C 


C 

c 

c 


c 


r 


SURPOUTINE  COMPKU*) 


THIS  SUHROOTINE  MANAGES  ALL  GUN  SYSTEM-ORIENTED  DATA  AND 
COMPUTATIONS,  the  INITIAL  ENTRY  POINT  READS  AND  PRINTS  SYST^ m 

data. 

common  X,Y.Z«XD0T. YDOT.ZnOT.XDD. YDD*ZOD.EGDOT.AGOOT,RSO. VP. 

. TF ,RNOS. VO.RMSO. THE AC.EDMAX. ADMAX .XMIN. XMAX. T I TLE (?0 ) .CALIR* 

.elmax.elgun 

COMmON/ENGPAR/P (9) 

COMMON/DYNtRR/DYNAM.nYNEM.nYNAV.nYNEV 

C0HM0N/HALIST/HARK.8ARK? 

COMMON/COVMaT/COV (?.?) 

COMMON/MEANS/  XMFANS  (f>)  .AMl  . am? 
real  PnS(6) tAVE(7) .C(9) 

equivalence  (x.posdn  t (x.c(i) ) 
integer  calir 

read  description  card.  CAlIPER.  GUN  RESIDUAL  DISPERSION  PARAMETERS 
» MUZZLE  velocity,  max.  TRACKING  RATES.  REACTION  Time. 
ROUNDS-PEH-RURST  AND  MAX  RANGE 

head  (?. 19.END=12)  TITLE.CALIH.SIGVO.SIGXD.SIGYD. VO.FOmaX. AOMaa. 

. treat,. RnOS.PMAX.ELMAX 
19  format  (?0A4/IS/4F10.0/(SF10.0) 

WRITF  <1.14)  TITLE.CALIR.VO.SIGVO.SIGXD.SIGYO.ELMaX.EOMAX. AOMflx. 

. TREAT. RMOS.RMAX 

14  format  < I0» .20a4/*0CALIHER  TYPE • . T94 . 15/ 

.•  MUZZLE  VELOCITY.  M/S • . TS4 . F8 . ?/ 

.*  STO.  OEV.  OF  MUZZLE  VELOCITY.  M/S • . T5G . FG . 4/ 

.*  STD.  DEV.  OF  X-COMP.  OF  RE9.  GUN  DISP..  R AD . • . TS7 . F9 . S/ 

.'  STO.  DEV.  OF  Y-TOMP.  OF  RES.  GUN  DISP..  RAD . • . T9 7 .F9 . M/ 

.•  Maximum  elevation,  rao. • .TSg.fp.a/ 

.•  maximum  elevation  RAIF.  R AD . /SEC. • . T56 . F6 . 2/ 

.»  MAXIMUM  AZIMUTH  RATE.  R AD . /SEC . • • T56 ♦ FG . 2/ 

.»  average  SYSTtM  REACTION  TIME.  SEC . • . T96 . F6 . 2/ 

..  ROUNDS  PER  RURST* .T5G.F4.0/ 

.«  MAXIMUM  EFFECTIVE  RANGE.  ME TERS • . TS4 . FG . 0 ) 

READ  IN  range  and  I IMF  OF  FLIGHT  USED  TO  COMPU 1 1 THF  HALLISTIT 
COEFFICIENT  KHAP 

read  (?.?0)  R.TF 
?0  format (?F10.0) 

RARK=VO/P-l .0/TF 
HARK?=RARK*BARK 

WRITE  OUT  all  Input  data. 

WRITFd.??)  R.TE.MARK 

??  FORwaTC  range,  m.  s'.F9.4.10X.i  TIME  OF  FllGHT,  SECS.  =*.FR.4. 
.lOX.*  PALLISTIC  COEFFICIENT. K-BAR  =«.FH.4) 

CALL  OYNGNl 
CALL  MFSIGl 

CALL  INITIALIZATION  PORTION  OF  VULNERABLE- AREA  SURROUTINE. 

CALL  VIILAPD(CALIR) 
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COMPXl 


OATF  s 7A10P 


1S/15/5S 


002A 

0025 

00?h 

0027 

002fl 

0020 


0030 

0031 

0032 

0033 
003* 

0035 

0036 

0037 
003R 
0030 
00*0 


00*1 


00*2 


00*3 

00** 

00*5 

00*6 

00*7 

004R 


00*9 

0050 

0051 

0052 

0053 
005* 

0055 

0056 


0057 


005H 


C COMPUTE  AUXILIARY  CONSTANTS  USED  FOM  BURST  KILL  PROBABILITY. 

RMS0=PMAX**2 

SIGV02=SI6VO**2 

SXRr)?sSIGXO**2 

SYG02=SIRY0**2 

RETURN 

ENTRY  COMPK(HPK) 

C THIS  F.NTRY  POINT  COMPUTES  BURST  KILL  PROBABILITY*  hpk*  GIVEN 

C THE  TARGET  POSITION. 

C THE  SINGULAR  POINT  AT.  THE  ORIGIN  IS  NOT  CONSIOEREU 
IF(X  .NE.  0.0  .OR.  Y .NE.  0.0)  GO  TO  100 
aGOOT  = 0.0 
EGDOT  = 0.0 
FLGUN  s 0.0 
GO  TO  1 
100  CONTINUE 

CALL  CTOP(C.P) 

C CALCULATE  THE  PARTIAL  OERIVATIVES 
CALL  6EUNC(A6i*EGl«2*f.l) 

ELGIJM  = EG  I 

VP=VO/ ( 1 .0*BARK*TF) **2  E 

SORSO  s SOHT(HSO) 

C MODIFY  ANY  MEANS  AND/OR  VAR.  WHICH  APE  EUNCTIO.S  OF  THE 
C ENGAGEMENT  PARAMETERS 

CALL  MESIG 

C perform  the  matrix  mutipl ication  necessary  TO 

C calculate  IHE  MEANS  AND  VARIANCES 

CALL  MATMUL 

C ADD  dynamic  pointing  error  to  The  MEANS  AND  VAR. 

C AND  CONVERT  TO  METER  ME ASURRMENTS  IN  THE  SLANT  PLANE 

CALL  DYNGUN 

AMI  s (AHI  ♦ OYNAM) *SORSO*COS (EL6UN) 

AM?  s (AM?  ♦ DYNEM)*SORSQ 

S3  = (C0V(1.1)*0YNAV)*RS0*IC0S(ELGUN) )**? 

S*  = (COV(2.?)  ♦ OYNEV)*RSO 
C THE  CORRELATION  COEFFICIENT 

RHO  = RSO* (COS(ELGUN) ) *COV< 1 *2) /SORT (S3*S4) 

C CALCULATE  THE  VAR.  DUE  TO  BaI  ISTIC  DISPERSION  ABOUT  TmR 

C FUTURE  POSITION  OF  IHE  TARGET 

XSAVE  = X 
X = X ♦ XD0T*TE 
CALL  CTOP(C,P) 

ADOT  = P(b) 

EOOT  = P|6) 

RVVsRSO/ ( VO*VP) 

SI  = < ADOT*RVV*COS<ELGUN) )**2*SIGV02  ♦ RSO*SXGD2 
S2=  ( EDOT*RVV)  **2*SIGVO?.RS')*SYG02 

C COMPUTE  INDUCED-CORRELATION  TERMS  - NOTE:  RHOB  IS  NOT  NOW  UREO  HUT 

C IS  INCLUDED  FOR  FUTURE  EXPANSION  OF  HPK. 

RHOB  s ADOT*EOOT*COS (ELRUN) *SIGV02*RVV**2/SQRT (S1*S?) 

C CALL  FOR  computation  OF  TARGET  VULNERABLE  AREA  GIVEN  TARGET 

C COORDINATES  IN  THE  ''POS'*  VECTOR. 

CALL  VULNARIPOS.VP. AVE) 
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oo«;<> 

c 

X « XSaVF 

CALL  FOR  COMPUfATION  OF  HURST  KILL  PROHABILITY.  HPk, 

c 

comrined  bias  AMO  standard  deviation  terms. 

00(^0 

CALL  BURPK  <AM1 »AM2.S1 . S? . S3 . S4 . RHOt HNOS t A VE ( 1 > »HPK) 

OOM 

RETURN 

006? 

1 

RPKsO.O 

0063 

return 

0064 

1? 

return  1 

006*; 

END 
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MAIN 


I'AIE 


7A10R 


Ib/lS/'-.S 


0001 

c 

c 

p 

SUBROUTINE  GFUNC(AG.ER»I«*) 

$ 

V, 

c 

c 

THIS  SUBROUTINE  CALCULATES  THF  PARTIAL  DERIVATIVES  OF  THF  GUN  POSITIU'’ 

c 

WITH  respect  to  THF  measured  OUANTITIES*  ANO  the  AZIMUTH  AND 

c 

elevation  PATES  OF  THE  RUN. 

Ono? 

C0MM0N/Z7TFLT/C0F (4> 

0003 

COMMON  OUMT(R) t FRO t AGO  * RSO » DUMY2 ( 1 ) . TFS»OUMY3» VO 

0004 

C0MMON/PARIS/F(?t6) 

OOO'^ 

COMMON/B aH ST/B ARK » BARK? 

OOOf) 

COMMON/ENGPAR/R* A .E*R0* AD»f  0*R00» AOOtEOD 

0007 

R2=P*R 

OnOR 

R02=RD*RD 

OnOR 

CE*COS(E) 

0010 

SE=STN(E) 

0011 

VT2=R2* ( aD*AO*CE*CE*ED*ED) ♦R02 

001? 

c 

VT*SORT (VT2) 

time  of  flight  FOUATION  linear  prediction  assumed 

c 

equation  assumes  the  projectile  FOLLOWS  THE  3/2  OPAR  LAW 

0013 

AAsPARK2*VT2 

r 

0014 

BB=?.0* {R*RD*BARk2*VT2*BAPK) 

r 

0015 

CC=R?*BARK2*4.0*P*RO«HARK*VT?-VO*VO 

r 

0016 

0D=?.0* (P2*B ARK *0*00 ) 

r 

0017 

EE=P? 

f' 

OOlR 

C0F(1)»PR/AA 

r 

OOIR 

COF (?) =CC/AA 

r 

0020 

COF (3) =OD/AA 

r 

0021 

COF (A) =FF/AA 

r 

002? 

CALL  TF1.T(  IF.MonO) 

0023 

TF?=TF*TF 

0024 

RF2=VT?*TF2*R2*?.0*R*RD*TF 

0025 

RF=50RT (RF2) 

00?6 

c 

IF(I  .FO.  3)  GO  TO  500 

calculate  the  terms  FOR  THF  PARTIAL  OERIVAIIVES 

0027 

raR*CF 

002n 

XO=Y*AO 

002R 

YO=bd*CE-P*EO*SE 

0030 

DEMON!  =Y  + YI)*TF 

0031 

XNUm1=XD*TF 

003? 

XLAMA=AT AN2 (-XNUM1 < DEMON 1 ) 

0033 

AG=A*XLAMA 

0034 

7=R*SE 

0035 

20=w*ED*CE»R0*SE 

0036 

xnum?=z*zo»tf 

0037 

T.J  = x^lUM?/RF 

* 

003R 

fG=aPSIN(T3) 

0030 

TFS  = TF 

0040 

RSO  = RF? 

0041 

T1=-XNUM1/0FM0N1 

004? 

T2=l .0/ ( 1 .0*T1*TI ) 

0043 

PAXn=-T?*TF/OE”ONl 
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GFUNC 


OATF. 


IA109 


li>/lS/S‘' 


00A4 

00*5 

OOAf. 

00*7 

on*M 

00*R 

0050 

0051 
005? 
0053 
005* 

0055 

0056 

0057 
005P 
0059 
00*0 
0061 
006? 
0063 
006* 

0065 

0066 
0067 
0069 

0069 

0070 

0071 
007? 
0073 
007* 

0075 

0076 

0077 
007fl 
0079 
OOPO 
OOPl 
OOP? 
OOP!? 
OOP* 
00R5 
OOP6 
OOP? 
OOPP 
00P9 
OORO 
0091 
009? 
0093 
009* 

0095 

0096 


PAYS  T?*T1/0EM0N1 
PAYnsPAY*TF 

PATFs-T?*XO*Y/ tnFMONl **?) 

PXOPsAO*CE 

PXOAOsY 

PX0F=-R*A0*SE 

PYRsCF 

PYORs-EO*SE 

PYOFn=-R*SE 

PYEsPYOEO 

PYDRO=PYR 

pYOfs-zn 

T*=1.0/SO«T(1.0-T3*T3) 

PE^sT*/RF 

PEZnsPEZ*TF 

PETFsPEZ*ZU 

PEHFs-PfZ«T3 

PZOEsYO 

pZORsFO»CE 

pzoRnsSE 

PZOEOsY 

PZRsPZDwn 

PZEsPZOEO 

PVTRsR* ( aO*AD*CF«CE*ED*EO) /vt 

PVTwOswn/V I 

PVTAnsR?*AO*CF*CF/VT 

PVTFs-«2*AO*A0*rE*SE/VT 

PVTFnsR2*E()/VT 

PRF VTsVT*TF?/WF 

PRFTFsl  VT?*TF*M*Pr))  /RF 

PRFRs (R*TF*RD) /PF 

PRFHnsTF*R/RF 

PAAVTs?,0*BARK?*VT 

PPBWs?.0*HO*RAPK? 

PBBWI)  = ?.0*R*HAM<? 

PHHVTs*.0*VT*PARK 

PCCRsPRBRD**.0*Rn*HARK 

PCCRns*.0*R*BARK 

PCCVT=?.n*VT 

POOMsPCCHD+?.0*RO 

PDORDs?.0*R 

PEERspnORO 

TF*sTF*** 

TF3sTF**3 

nEMON3  = *.0*AA*rF3*3.0*HR*TF?  + ?.0*Cr.*TF  + nn 

PIFRs- (TF3*PRRR*TF?*PCCR*TF*P00R*PFER) /OEMON  I 

PTFRO=- (TF3*PBHRn*TF?*PCCRO*TF*POORD)/nEMON3 

PTFVT  = - (TF**PAAVT*IF3*PRBVT»TF?*PCrVT)  /OF.MOi''  i 

PTFRsPTFR*PTF  V I *PVIR 

PTFRnsPTFRU*PTFVT*PVTRn 

PTFADsPTF  VT*PVT  All 

PTFEO=PTFVT*PVTEn 

PTFE=PTFVT*PVTE 
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15/lS/SS 


0097 

009fl 

0099 

0100 
0101 

0102 

0103 

0104 
OlOG 
010(S 
0107 
OlOO 

0109 

0110 
0111 
0112 

0113 

0114 
01 IG 
01  Ifi 
0117 
OtlG 

0119 

0120 
0121 
0122 
0123 


PRER=PRFR*PRFVr*PVTR»PRFTF*PTFR 

PHFRD«PRFRO*PRFVT*PVTRO*PHFTF*PTFHn 

prfao*prfvi*pvtad*prftf*ptfad 

PRFFn=PRFVT*PVTEO*PRFTF*PTFEO 

PRFr=PRFVT*PVTE*PRFTF*PTFF 

: CALCULATE  THE  PARTIAL  flEPIVATIVES  USING  THE  CHAIN  RULE. 

F(l»l)  = PAXO*PXnR*PATF*PTFR*PAY*PYR*PAYO*PYDR 
Fdt?)  = 1.0 

F(1.3)  = paxo*pxoe*patf*ptfe+pay*pye*payo*pyoe 

F(1.4)  = PATF*PTFRO  ♦ PaYO*PYDRD 
F(l.G)  = PAXO*PXDAO  ♦ PaTF*PTFAO 
F(l.h)  * PATF*PTFEU  ♦ PAYO*PYOEn 

F(2.1)  * PEZD*P70R  ♦ PETF*PTFR  ♦ PFRF*PRFR  ♦ PEZ*PZP 
F(2,2)  = 0.0 

F(2.3)  = PEZD*PZOE*PETF*PTFF*PERF*PHFE*PEZ*PZF 
F(2.4)  = PEZO*PZnRO*PETF*PlFRO*PERF*PRFRI3 
F(2,5)  = PETF*PTFAO  ♦ PERF*PRFA0 
F(2.6)  = PEZn*PZnEO  ♦ PETF*PTFED  ♦ PERF*PRFEG 
: CALCULATE  THE  RATES  OF  THE  GUN  IN  AZIMUTH  AND  ELEV. 

DTFDT  = PTFR*RD»PTFR0*R0n»PTFA0*A0D+PTFE0*E0n*PTFE*FO 
A6D=AD-(OEMON1*XO*OTFOT-XNUM1*(YD*YD*OTFOT) )/(DEMONl »*2*XNUMl**2) 
ORFDT  = PRFR*RD*PRFRO*RDn*PRFAD*AOO*PRFEO*EOO*PRFE*Fn 
E60=T4* (RF*ZD* ( I . ♦OTFOT ) -XNUM2#0RF0T  ) /RF/RF 
RETURN 

: NO  FIRE  CONTROL  SOLUTION 

1000  RSO  = 1.0E12 
RETURNl 

500  RSQ  * RF? 


RETURN! 

END 
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OoOl 


000? 

0003 

0004 

0005 
OOOG 
0007 
OoOft 

0009 

0010 
0011 
0012 

0013 

0014 

0015 

0016 
0017 


C 

C 

C 

C 

C 

C 


SUBROUTINE  CTOP  (CfP) 


COMPUTE  THE  POLAR  COORDINATES  FOR  A POINT 
LOCATED  BY  CARTESIAN  COORDINATES  (C(D) 
dimension  C(9) tPCR) 

SS  r C(1)*CI1)  ♦ C(2)«C(?) 

P<1)  = SORT(SS  ♦ C(3)<»C(3)) 

S * SOPT(SS) 

SDOT* (C ( 1 ) ^C (4) ♦C (2) ♦C (S) ) /S 

P(2)  = ATAN?(C(2)  •C(l)  ) o 

IF(P(?)  »LT»  0*0  .AND*  C(l)  .LT.O.O)  P(?)  = P(2)  ♦ 6.2B31H 

P(3)  = ATAN?(C(3) *S) 

P(4)  = (C(1)*C(4)  ♦ C(2)<>C(S)  ♦ C(3)*C(6)  )/P(l) 

P(5)=(  C(l)*C(b)-C(4)«C(?) )/SS 
P(6)  = (C(6)  -C(3)*P(4)/P(1) )/5 

P(7)  = (C(4)*C(4)*C(5)*C(S)^C(6)‘*>C(6)^C(1)*C(7)^C(?)<^C(B)*C(3)*C(9) 

1 -P  ( A)  *P  (4)  ) /P  ( 1 ) 

P(H)  s (C(B)*C(1)  - C(?)*C(7)  -?.0«S*SOOT«P(5) )/SS 

P(9)=(C(9)-(P(1)*(P(7)*C(3)^PC4)*C(6))-P(4)*P(4)*C(3))/ 

I (P(l)^P(l) )-P(6)*SD0T)/S 


RETURN 

END 


non 
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DATE 
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Ib/lS/HS 


0001 


ono? 


0n03 

0004 

ooos 

oooA 

0007 

OOOfl 

0009 


0010 

0011 

0012 

0013 

0014 

oois 

0016 

0017 


001  s 
OOlR 


0020 

0021 

0022 


SUHROUTINE  plot (RMIN, BMAX«B. LABEL) 

C 

C THIS  SUBROUTINE  PLOTS  THE  PK  CONTOURS  AND  SYSTEM  CONSTRAINTS  IN 

C THE  PLANE  WITH  X AS  THE  HORIZONTAL  AXIS  AND  B AS  THE  VERTICAL 

C AXIS.  IF  BMJn  UMAX  AND  B HAVE  THE  VALUES  OF  THE  RESPECTIVE 

C VARIAHLES  of  Y*  THE  PLOT  WILL  BF  A "GROUND  PLOT".  IF  THE  DUMMIES 

C BMlNt  BMAX  and  R CONTAIN  THE  ASSOCIATED  VALUES  FOR  Z.  THE  PLOT 

C WILL  HE  IN  THE  VERTICAL  PLANE  ("SKY  PLOT"). 

C 

C THE  PROGRAM  WORKS  BY  INCRFMENTING  THE  DUMMY  VARIABLE  B HETWFEN 

C HMIN  and  BMaX  during  THE  PROCESS  OF  SWEEPING  THF  PLOT.  THE 

c variablf  in  the  calling  list  for  B will  be  changed  accordingly. 

C THUS  allowing  for  the  freedom  of  choice  to  either  SlvEEP  Y OR  2 

C MERELY  BY  PLACING  THE  VARIABLE  NAME  IN  THE  LISl. 

C 

COMMON  X.Y.Z.XDOT.YOOTtZPOT.XOD.YOD.ZOD.FGOOT. AGOOT.RSO. VP. 

. TE.RNnS,VO,RMSO.TREAC»EDMAX»ADMAX.XMIN.XMAX,TITLE(20) .CALin. 
.FLMax.ELGUN 

common  /SPLOT/  pk.char.mchar.npk 
C0MM0N/EMGPAR/P(9) 

eouivalence  (C(i).x) 
dimension  C(9) 

real  AAA  ( IS)  .PSAVEdOl  .SI ) fPK  (10)  .LABEL  (4)  ,H  ( 1 ) 

INTEGER  *2  CHAH ( I 0 ) . MCH AR ( 4 ) . IPLOT ( 1 0 1 . 5 1 ) . BL ANK/ • •/ 

INTEGER  CALIB 
C 

C PERFORM  INITIALIZATION  FOR  THE  PLOTTING  SUBROUTINE  "SPLOT". 

C 

AAA ( 1 ) =XMIN 
AAA  (2) =XMAX 
AAA (3) =BMIN 
AAA (4) sHMAX 
AAA (S) s-S.O 
AAA (6) =-3.0 
CALL  SPLOTO(IPLOT) 

CALL  SPLOTI(AAA) 

C 

C GET  THE  increments  FOR  THE  AXES  FROM  SPLOT. 

C 

OX=aaA (7) 

0B=AAA  («0 
C 

C store  zeros  into  a matrix  used  to  STORE  PK  AT  EACH  POINT  ON  GRID. 
C 

DO  12  IX=1.101 
00  12  18*1*51 
12  PSAVE(IX.IB)=0.0 
C 

C initialize  X EOP  SWEEP  FROM  LEFT  TO  RIGHT. 

C 
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IS/IS/SR 


0023 

0024 


002*i 

002#. 


0027 

0020 


002P 

0030 


0031 


0032 

0033 

0034 


003R 

003#. 

0037 

003H 

0030 

0040 

0041 


0042 

0043 


0044 

0045 

0045 

0047 

004« 


0>  1 IX=ltlOl 
C 

C initialize  a "PART"  VALUE  FOR  PK. 

C 

R( llsHMIN-On 
CALL  COMPK(PKS) 

C 

c initialize  the  0 variable  for  sweep  from  bottom  to  top. 

c 

BID  =BMIN 

no  2 iHsi.si 
c 

C COMPUTE  PK  AT  THIS  POINT  AND  SAVE. 

C 

CALL  COMPK(PKC) 

PSAVE(IX,IH)=PKC 

C 

C TEST  IF  TARGET  IS  OUTSIOE  MaXIMIUM  EFFECTIVE  GUN  RANGE. 

C 

IF  (RSO.LE.RMSO)  GO  TO  3 
C 

C store  an  "R"  if  first  TIMF  max.  range  excefueo. 

c 

IF  (IM.EO.l)  GO  TO  1 
IPLOT ( IX. IH) sMCHAR(4) 

GO  TO  1 

r 

C COMPUTF  RANGE  AT  WHICH  TARGET  WOULO  ME  IF  IT  CONTINUED  ON  ITS 

C.  PRESENT  COURSE  FOR  A TIMF  EOUAL  TO  THE  AVERAGE  SYSTEM  REACTIOrj 

C TIMF.  TRFAC. 

C 

3 XS  = X 

X = X ♦ XOOT»TREAC 

IF  ( X .FQ.  0.0  .AND.  r .FU.  0.0  ) GO  TO  1000 
CALL  CTOP(C.P) 

CALL  RFUNC(  AOUM,EDUM.3.E-10  00) 

1000  X a XS 

RNSO  a RSO 
C 

C TEST  IF  NEW  RANRF  IS  HEYONO  MAX,  SYSTEM  RANGE  AND  PLOT  A "T". 

C 

IF  (RNSQ.LE.RMSO)  GO  TO  4 
lPLOT(IX,IH)rMCHAR(3) 

C 

C test  if  either  max.  AZIMUTH  RATE  OR  max.  ELEVATION  OR 

C elevation  hate  are  EXCEEDED 

C 

4 IF  (AHS(AGOOT) .LE.AOMAX)  GO  TO  5 
IPLOT ( IX. IB) aWCHAR (2) 

GO  TO  G 

S IF(AHS(FLGUN)  .GF.  ELMAX)  go  to  2000 
IF  (APS(FGOOT) .LF.FOMaX)  go  to  5 
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Or>49 

OOSO 


On«?l 


0092 

0093 

0094 


0095 

0096 

0097 


OOBfl 


0099 


0060 


0061 

006? 

0063 

0064 
0069 
0066 

0067 

0060 

0069 

0070 

0071 
007? 

0073 

0074 
0079 


2000  IPL0T(IX»I8)smChaR(1) 

60  TO  9 
C 

C SET  SIGN  OF  PK  GRADIENT  TO  SS  FOP  TEST. 

C 

6 SS=SIGN( 1 .OtPKC-PKS) 

C 

c SCAN  THF  PK  VECTOR  TO  SFE  IF  PKC  H«S  CROSSED  ONE. 

C 

DO  7 IP=1,NPK 
PKDePK (IP) 

IF  (.NOT, (SSO(PKD-PKS) .GF.0.0  .AND.  SS* (PKO-PKC ) . LE . 0 . 0 ) ) GO  TO  7 
C 

C IF  PKC-PKC  straddles  PK(IP).  plot  THE  APPROPRIATE  CHARACTER. 

C 

iplot(Ix,ib)=char(ip) 

GO  TO  « 

7 CONTINUE 
C 

C INCRFMFNT  B UP  RY  ONE  LINE. 

C 

fl  R(1)=P(1)*0R 
C 

C MAKF  PKS  EOUAL  TO  PKC  FOR  NEXT  PASS. 

C 

? PKSsPKC 
C 

C increment  X RY  one  to  THF  WIGHT. 

C 

1 X=X»DX 
C 

C THE  SCANNING  OF  THE  GRID  FROM  HOTTOM  TO  TOP  ANO  LEFT  TO  RIGwT  Is 

C complete.  NOW  SCAN  FROM  LEFT  TO  RIGHT  AND  HOTTOM  TO  TOP  TO  fIlL 

c IN  THF  lines  as  needed. 

C 

DO  9 IH=1,S1 
PKSsO.O 

DO  10  IXrl.lOl 
PKC=PSAVF (IX, IR) 

IF  (PKC.F(J.0.0  .OH.  PKS. EO. 0.0)  GO  TO  10 

IF  ( .NOT. ( IPLOT ( IX, IB) .EO.HLANK  .OR.  IPLOT ( I X , I B ) . EO. MCH aR ( 3 ) ) ) 

. 60  TO  10 

S=SIGN ( 1 ,0 ,PKC-PKS) 

DO  11  IP=1,NPK 
PKOsPK ( IP) 

IF  ( .NOT. (S* (PKD-PKS) .GE.0.0  .AND.  S* (PKO-PKC ). LE . 0 . 0 ) ) GO  TO  11 
IPLOT ( IX, I») =CHaR (IP) 

GO  TO  10 
11  CONTINUE 
10  PKSsPKC 
9 CONTTNUF 
C 

C PRINT  OUT  THE  ''FOOTPRINT" 
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C 

0o7ft  CALL  9PL0T3( IPLOT*AAA*20HISO-PK  CONTOUR  HAP 

. 16H00WN  RANGE*  H,  ,1) 

0077  RETURN 

007R  ENO 


tLABEL. 
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onoi 


ooo? 

0003 

0004 

0005 

0006 

0007 


OOOB 


OOOQ 

0010 

0011 

0012 

0011 


0014 

0015 

0016 
0017 


C 

C 

C 

C 

r 

c 

c 


c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


=rss^sss=ssss=s3xss=sssss=s==;sssssss==ss==sss3sxs=ssssss 

SUBROUTINE  SPLOTO ( I PLOT ) 

33 SS3SS3.3SS 3 33 33 Ss 33 3 3333333333333333333333333X333333333 

SPLOTO  initializes  THE  'PICTURE  MATRIX'  HY  CLEARING  IT  TO  BLANKS.  P0000020 
THE  PICTURE  MATRIX  HAS  lOl  COLUMNS  CORRESPONDING  TO  THE  101  OOOOOOIQ 

HORIZONTAL  HAMMER  POSITIONS  ON  THE  PRINTER.  THE  PICTURE  hAToiX  HAS00000040 
51  PnwS  CORRESPONDING  TO  THE  VERTICAL  SPACING  ON  THE  PRINTER.  OOoOOOSO 


integer  *?  IPLOTdOl.51)  tMLANK/'  «/ 
DO  1 TYxI.'dI 
00  I TXxl.lol 
1 IPL0T(TX.IY)3HLANK 

return 


0000006(1 
O0O00070 
onoooOHO 
01)000000 
M)0001 00 


entry  SPLOTKAAA)  OOOOOllO 

SPLOTl  INITIALIZES  THE  WORKING  PORTION  OE  THE  USER'S  AAA-VEfTOW.  (>0000120 


real  AAA(IS) .EPSLON/5.0F-4/ 


O0000140 


USER  MUST  PROVIDE  A VECTOR  OE  LENGTH  15  *4  WORDS. 

USER  provides  PARAMETERS  IN  EIRST  6 WORDS.  THE  LAST  O ARE 
USED  HY  THE  PROGRAM  EOP  STORAGE  OE  INTERMEDIATE  RESULTS. 
AAA(1)bXMIN 
AAA (2) xXHAX 
AAA(3)aYMlN 
AAA(4)xYHAX 
AAA(S)xTXRWR 
AAA(6)sTYPWR 


00000150 

00000160 

00000170 


AAA(7)c(AAA(2)-AAA(1)) /lOO. 
AAA (A) 3 (AAA (4) -AAA (3) ) /SO. 
AAA ( 15) 30. 

RETURN 


00000240 

00000250 

M0000260 

(.0000270 


entry  SPLOT?(IPLOT.CHAR.AAA.XV.YV)  00000320 

SPLOT2  determines  THE  PROPh R POSITION  IN  THE  GIVEN  PICURE 
MATRIX  OE  IHF  PAIR  (XV.YV)t  USING  THE  SCALING  INFORMATION 
given  in  THE  AAA-VKCTOR. 

IPLOT  is  THE  )*ICTUHE  MATRIX  1HE  (XV. YV)  POINT  IS  TO  ME  PLOl IM)  ON.OnnOOlin 
CHAR  IS  THE  PLOTTING  CHARACTER  TO  ME  USED.  (I NT E(iERa:> ) 

AAA  TS  The  PLOTTING-PARAMf TFR  VECTOR  FOR  THE  FUNCTION  mEING 
PLOTTED  (SEE  AMOVE). 

XV  and  YV  APE  the  COORDINATES  OF  THE  POINT  TO  ME  PI  OTTED  MY  SPLOT200000170 


integer  *2  char 

INCREMENT  A POINT  COUNTER. 
NPsAAAdS) 

NPsNP* 1 
AAA (IS) 3NP 


00000300 

O0O00410 

00000420 

00000430 


4 
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IS/ 1 s/ss 

c 

r.ei  THf  DELTA-X  ANO  DFITA-Y  FOR  THIS  FUNCTI0>4. 

OOIR 

HXsaaA  (7) 

OO00045O 

onl9 

nY=AAA (H) 

110000460 

c 

lEST  IF  PAST  IMTTIAL  START-UP  PFRIOD  OF  1 POINTS. 

oo?o 

IF  (MP.GT.3)  go  TO  61 

(i0niU)4Mn 

c 

SAVF  THE  NEW  pniMTS. 

Oft?l 

AAA (NP«8) sAV 

oononson 

00?2 

AAA (NP^ 1 1 ) sYV 

noooosin 

c 

IF  WF  NOW  HAVE  1 POINTS#  Si  APT  IMF  PLOT. 

•> 

oo?i 

IF  (NP.FO.J)  GO  TO  60 

onoonsTO 

00?4 

return 

onoons40 

c 

SAVE  NEW  POINTS  AND  PUSH  (illl  OJ  I>FSr  (TNF  . 

0O?S 

51 

1 1 *<? 

00?^ 

00?7 

AAA (9) 3AAA ( lU) 

00?fl 

AAA ( 1 0 ) xA AA ( ] n 

60000570 

0029 

AAAdl  )=XV 

4H1O00S8O 

0030 

AAA ( 1 2) sAAA (13) 

•*0000590 

0031 

AAA(13)=AAA( 14) 

00000600 

0032 

AAA (14) xYV 

00000610 

0033 

GO  TO  63 

00000640 

0034 

60 

Ll*l 

003S 

L2=? 

c 

COMPUTE  A DIRFCTION-SFNSITIVF  UELTA-X. 

003f> 

61 

OYTxARS ( AAA (L2^ 1 1 ) -AAA (1  1 ♦ 1 1 ) ) 

0037 

CORR*OX*nY/AMAXl (OY#nYT) 

003R 

COPRxAMAXI (CORR#0.01*nX) 

0039 

OIPECT  = STGN(COPPt  AAA(L2*H)-AAA  (Ld8)  ) 

c 

initialize  an  x-scanner  to  TRAVFL  between  points  L1  and  L2. 

ornOOYOO 

0040 

XSCANsAAA(Ld8) 

0041 

IF  (L1.E0.2)  XSrAN=XSCAN40lRECT 

c 

rOMPUTF  ANO  test  values  USFO  HY  newton  interpolation  FORMULA. 

flOOOO  7H0 

004? 

nENlxAAA (10) -AAA (9) 

0OO00790 

0041 

0EN?»AAA (11) -AAA ( lU) 

0 n 0 0 0 »( 0 •) 

0044 

nENl=AAA(ll)-AAA(9) 

OOO0081 0 

0045 

IF  ( ABS(hENl) ,LT.EPSL0N#APS( AAA (10) ) ) GO  Id  64 

oonnortpo 

004(S 

01 « (AAA (13) -AAA (12) ) /OFNl 

OOOOOH  10 

0047 

GO  TO  65 

noo00840 

004R 

64 

01*0. 

onnoOHSO 

0049 

65 

IF  ( ABS(DEN?) .L r.EPSl ON«ABS(AAA( 10) ) ) GO  TO  66 

00000860 

0050 

02«  (AAA(14)-AAAn3))  /DEN? 

00000870 

0051 

GO  TO  67 

OOOOOHHO 

005? 

66 

02  = 0. 

00000890 

0053 

67 

IF  (AHS(nENl) ,LT.FPSLON*AMS(AAA(10) ) ) (U)  TO  6H 

on0O0900 

0054 

noix (D?-ni ) /OENi 

000009I 0 

0055 

GO  Tn  69 

on O00920 

005R 

68 

DDlxO. 

00000910 

c 

USE  NFwTON'S  formula  TO  COMPUTE  A Y-VALUF  CORRE SPONO I NG  TO 

XSCAN. 

00000940 

0057 

69 

iXxi .5* (XSCAN-Aaa ( 1 ) ) /ox 

c 

TEST  IF  TT  IS  OFF  SCALE. 

005R 

IE  (tx.lt.1  .or.  ix.GT.inu  return 

nooo0770 

0059 

YSCANxAAA  ( 1?)  ♦(XSC4N-AAA(9))«(I)1*  (XSCAN-AAA  ( 10)  ) ♦OD)  ) 

r 

COMPUTE  rORRESPONOiNG  VERTITAL  SPACE  POSITION  ON  PRINTER. 
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15/lS/SK 


00^0 

0061 

0062 

0063 

0064 

0065 

0066 


0067 

0068 

0069 

0070 


I 


0071 

007? 

0073 

0074 
007*; 
007f» 

0077 

0078 

0079 
OOPO 
0081 
008? 

0083 

0084 
0086 
0086 

0087 

0088 

0089 

0090 

0091 
009? 

0093 

0094 


C 

C 

C 


C 

C 

C 


C 

C 

C 

C 

c 

c 

c 

c 

c 

c 

c 

c 


I Y=1 .64 (YSCAN-AAA (3) ) /DY 
TEST  IF  OTf  SCALF  A^»0  CORPfCT. 

IF  (lY.LT.U  lYsl 
IF  (TY.GT.51)  lYr5J 

STORF  A plotting  CMARaCTFR  IN  ThF  DFSIREH  POSITION. 

IPLOT ( IX# I Y) cCHAR 

INCRF'^rNT  X-6CA^♦N^P  AMD  TFM  IF  PAST  POIMf  L?. 
ASCAN*XSrAN4|)IRF  f I 

IF  (nTMFCT*(XSCAN-AAA(L?^H)  ) ,Lf  .0.0)  GO  10  f>9 
IF  PAST  POINT  AND  NOT  IN  INITIAI  PtDIOU#  oMURN. 

IF  (LH.EO.J)  RtTURN 

JUST  PLOTTED  PLTWEEN  THF  FIRST  ? POINTS#  NOW  PLOT  PFTWFtN  THE 
?-ND  AND  1-RI)  POINTS. 

Ll=? 

L2=3 

GO  TO  63 


00000970 

00000990 
00  0.0 100  0 

fR)O010?O 

oooolo  ro 

<H)00  1 OMO 
MOO  ()  1 090 
01)001 1 on 
oonoilio 
00001120 

OOOOl 140 

noooiiso 

OOOOl 1 70 


FNTPY  SPLOT3C IPLOTf AAA#TITLF  fVLARfHLAHfO) 


M0001210 


SPLOTT  PRINTS  OUT  THF  PTCTHWE  MATRIX  USING  THI  LAHFLS  AND  OooOl l90 

THf  scaling  information  GIVhN. 

TITLF  is  a 20-CHAHACTEW  HOLLFWITH  FIFU)  GIVFN  hy  THF  USER  TO  L AHf L onoO 1 ??0 
THE  TOP  OF  THE  PICTURF 

VLAH  is  a 16-CHARACTFP  HOLLFRITH  field  GIVFN  HY  THE  USER  TO  LAHFL  000l»l?40 
THF  vertical  scale. 

HLAH  TS  A 16-CHARACTFR  HOLLFRITH  FIFU)  GIVEN  RY  THE  HSEP  TO  LaMFL  0on01?60 
THE  HORIZONTAL  SCALE. 

0 IS  THE  OUTPUT  DATA  SET. 


peal  TITLE(5)#VLaR(4) fHLAHCA) »XLAH(ll) 
INTFGFR  0 
WRITF  (0#?)  TITLF 
? format  (tlt,50XSA4) 

HIGOY*10.0»*1F1X ( AAA (6) ) 

YLAR.AAA (4) •HIGDY 
RlOOYaAAA («) OMlGnY 
IP=0 

no  ♦ j«i#5i 
I Y«s?-J 
1P«1P^  J 

IF  (J.NF.23)  go  to  ?6 

WRITF  (n#21)  VLAB# ( IPLOT ( I X • J Y)  #lX  = l#lOn 
GO  TO  31 

?6  IF  (J.NF.24)  GO  TO  27 

l=aaa  <«>) 

wRiTF  (0#22)  L# (IPL0T(IX#1Y) #lXrl,lOl) 

GO  TO  31 
?7  IF  (lP)St6#6 

6 TP=-S 

WRITF  (0.7)  YLAR. ( I PLOT  ( T X 1 1 Y ) # I X = 1 # 1 0 1) 

7 format  (*  ••F16.?#^H  *101AI) 

GO  TO  31 

5 WRITF  (0#fl)  (IPLOT(IXflY) tIX=l#ini) 


00n0l?90 
OOO01300 
0 1)001310 

nnoOl 330 
nnnoi 340 
on 00 1 3S0 

TiOMOI 
M0001370 
(iMOOl  380 
00001390 
n n 0 0 1 4 0 0 
MOOOl 410 

M r)ooi4?o 
^0001430 
00001440 
O0O014S0 
00001460 
00001470 
000014HO 
r,0n01  490 

0000  ism 

nnoo 1S?0 
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is/ls/^^ 


009S 

8 

FOMMAT  ( • • f 1 7X • 1 • f 10) Al) 

00<)6 

11 

IF  (MOIM  Jf 6)  ,F0. 1)  (iO  TO 

00P7 

WMITF  (0.?9) 

OOPfl 

?9 

FOMMaT  ( • ♦ • . llHKf  • I • ) 

009P 

GO  TO  4 

0100 

?« 

WWITF  (0.30) 

OlOl 

^0 

F0HHAT(»4**19X.  ?0(Sh  ♦)) 

010? 

IF(  J.FQ.  1 .OM.  J .F(i.  SI)  WMlTFCOf' 

0103 

50 

FOPmaT ( • ♦ • • IftXf *-• .20 (5h — --  )) 

0104 

4 

VLAH»YL4H-dIG0Y 

0105 

Rl6OXxl0.0<»*IFTX(AAA(5)  ) 

0106 

XLAR(l)sAAA(l)*R16DX 

0107 

TENOXsIO.^Aaa (7)*BIG0X 

0108 

DO  9 J=lflO 

OlOQ 

9 

XLAR(j4l)»XLAR(.J)  ♦TENDX 

0110 

L*Aa A (5) 

0111 

WRITE  (0.10)  XLA^fHLAR.L 

01  1? 

10 

format (•  •flRXf»|».lO(*....*....IM/ 

1 14H  >iULT.  PV  10*»I2) 

0113 

PETUWM 

0114 

PI 

format  (•  * .4A4 , 1 X » 1 • f 1 0 1 A 1 ) 

0115 

format  (t  I.14M  mult,  my  I0*»l?fix»l 

0116 

PNC) 

• • • I IF 10*?» 


• « lOlAl) 


l‘«M)U  1 S U 
onon  1 

onoOl S U 
iinno  1 s Pi 


oonoiSAO 
onooisso 
oonoisfio 
oonOlSTo 
oonoiSRO 
oonoiS90 
00001600 
0 0 n 0 1 M n 

t 

0000)610 
0 00  0 16<»0 


67 


FOOTRAM  TV  G LFVFL  21 


MAIN 


DATE 


74109 


lb/lS/S5 


c 

p 

Onoi 

subroutine  MATMUL 

C 

c 

c 

This  subroutine  performs  the  matrix  MULTIPLICaTIOS  necessary 

c 

TO  Calculate  the  means  and  var,  of  the  gun  error 

oon? 

COMMON/COVMAT/  COV(2t2) 

0003 

COMMON/PaRTS/  F(?.6) 

0004 

common/  sigma/  SIG(6t6) 

OOOG 

COMMON/MFANS/  XM(6)«AM(2) 

0006 

DIMENSION  TEMP(2»6) 

0007 

DO  10  I=lt2 

OOOft 

AM(I ) =0. 

0009 

DO  10  J=lt6 

c 

CALCULATE  THE  MEANS 

0010 

AM(I)  = aM(1)  ♦ F(ItJ)*XM<J) 

001  1 

SUM  = 0.0 

ool? 

DO  S K=l,6 

0013 

5 

SUM  = SUM  ♦ F(I.K)*SIG(Kt J) 

0014 

10 

TEMP(I.J)  = SUM 

OOlS 

no  20  I=l»2 

0016 

DO  20  J=T»2 

0017 

SUM  = (1.0 

0018 

no  IS  K=lt6 

0019 

15 

SUM  s SUM  ♦ TEMP(I.K)*F(J»K) 

c 

calculate  the  VAR. 

oo?o 

COV(I.J)  = SUM 

0021 

?0 

COV(J.I)  = SUM 

002? 

RETURN 

00?3 

END 
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main 


DATE 


M109 


1 


t 


onoi 


Oon? 

0003 


0004 

onos 


onoA 

0O07 

Ooon 

0009 

0010 

0011 


ool? 

001  3 


0014 

0015 


• 0016 
0017 


c 

subroutine  MESIBl 

C 

COMMON/MFANS/  XM  ( ft ) t AMOlIM  ( 2 ) 
common/sigma/  SIR(ft.ft) 

C head  THE  GUN  SE^'SOH  ERROR  OATA 

c The  order  of  the  varihales  in  the  sensor  mean  mEaSURFmFNT 

C ERROR  VECTOR  (XM(D)  IS  : 

C * RANGE 

C * A7IMUTH 

C * El EVaTION 

C * RANGE  HATE 

C * AZIM(JTH  HATF 

c * elevation  Rate 

MS? 

REAOIM.IOOO)  (XM(1) .Isl.ft) 

c The  input  sig(I*j)  array  contains  the  stanoaho 

c deviations  in  the  diagonal  te‘ ms  and  the  correlation 
c coefficients  in  the  off  diagonal  terms,  the  order  of  the  hors  of 

C THE  matrix  (and  the  cols.)  CORRESPOND  TO  THE  KOl.LOWiNG  OHOER  OF  THF 

c Engagement  parameters: 

C 1 sr>  RANGE 

C ? S=>  azimuth 

C ? s=>  elevation 

C 4 ==>  range  rate 

c s ==>  az.  rate 

c ft  ==>  FLEW.  RATE 

c The  sig  Matrix  rul  be  a symetricai.  matrix  and  it 

C IS  REAO  ONt  ROW  at  a TIMF.  THE  INPUT  DATA  SHOULD  BE  SET 

C UP  WITH  ONE  ROW  OF  THF  MATRIX  ON  A CARO. 

RE AO (M, 10 00)  ((SIG(J.I).Isl,ft),J=l.ft) 

1000  FORMAT (ftFl 0.0) 

WRITE  (1.100) 

100  F0RMAT(*  SENSOR  PRRORS*/*  MEAN  MFASURENT  ERRORS  XM(t)»/) 

WHITE  (1.200)  (XM(I) .Isi.ft) 

?np  format  (ftX  » •RANCiE.M'  .T30.FB.S/ftX.  • A/ 1 MUTH  i RAD  * . 1 30.FB.S/ 

.ftX. « elevation. RAO* . T30.FH.S/ftX. * RANGE  RATE. M/S* . T30.FH.S/ 

.ftX.'AZ.  PATEtRAiVS* .T30.FH.S/ftX, *EL.  RATE .RAD/S* .T30.EB. 5/) 

WRITE  (1.300) 

300  FORMAT!*  matrix  of  STaNOARU  DEVIATIONS  AND  CORRELATION  COEEFlCirwr 

.S'/ 

.*  THF  UNITS  ARE  the  SAMF  AS  FOR  THE  MEANS*// 

. IftX. 'RANGE  * .ftX . * azimuth* .7X.*ELEV.*.7X.*R.  RATf'.ftX.'AZ.  RAIL*. 
.SX.'EL.  RATE*/  ) 

WRITE  (1.400)  ( (SIG( J. I ) . 1=1 .ft) . Jsl ,6) 

400  FORMAT!*  RANGE • .ftX.FR.2.ftX.FH.3.4 (SX.EP.3) /*  AZIMDTH  • . ft ( 4 X , FR . a ) / 
.*  Elevation* .3X.F9. 4. 5(4x.F9, 4)/*  r.  rate* .ft(sx.EB.3)/ 

.*  AZ.  RATE  * .4X.F9.4.5 (4X.F9.4) /*  EL.  RATE*. 4 .E9 . 4 . S ( 4 X . E9 . 4 ) ) 

C Change  the  sto.  oev.  into  var. 

00  5 T=1.S 
K = T ♦ 1 
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FOpTRAN  IV  6 LEVEL  21 


MFSIfil  DATE  = 7A109 


OOlfi 

no  S J =K«9 

0019 

SIGd.J)  = SlG(I.J)*Sir,(I,I)*SlR(J.J) 

On?o 

siGcj.n  * siG(i.j) 

on?i 

s 

CONTINUE 

OOP? 

00  b I=lt6 

00?3 

6 

1>IG(I,I)  = SIGdfl)**? 

002A 

RETURN 

On?S 

END 

1S/15/SS 
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main 


DATF  = 74ia<) 


1 S/IS/S^ 


FORTRAN  IV  0 LEVEL  21 


0001 


000? 


0003 

0004 


0005 

0006 
0007 
0000 


Ooog 

0010 

0011 

0012 

0013 


C 

C 

C 

C 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


c 

c 

c 


c 

c 

c 

c 

c 

c 

r 

c 


sxrsssr=sssssz7ss=ssxrsss£?s?rsrsssssss=sssrsBxssssxsss= 

subroutine  VULAPD<CALIB) 

SXSX3SS3SSSSSSSSSSSSSSSSSsSSSZSSSZ=SSSSSSSSSSSSSSSSSSsS= 


THIS  SUHROUTINF.  COMPUTE^  TARGET  VULNERABLE  AREA  FOR  PO I NT -CONT ACT 
Projectiles. 

vulnerable  area  is  a FUNCTION  OF  TARGET  TYPE*  bHELL  rALlHFR, 
STRIKING  VELOCITY  AND  STRIKING  DTRFCTION. 

THIS  ENTRY  POINT  SERVES  TO  SELECT  ThF  OaTA  FOR  THE  SHELL  CAl IHPR 
BY  "CALIR**. 


implicit  INTEGER  (E) 


nonofipon 


CALlBxl : 

CALIH*?: 

CAlIB.31 

CALIBsa: 

CALIHzS: 

CALIHz6: 

CALIHc7: 

CAlIBzO: 

CAlIBxR: 

CALlBalOl 

CALlBzll ; 

CALlBzl?: 

CALIHz13: 

CALIHzIA! 

CAlIB*15: 


20MM  HF  PROJ.*  KK  KILL. 

30MM  HF  PROJ.*  KK  KILL. 

40mm  he  PROJ,*  KK  KILL. 

57MM  he  PROJ,*  KK  KILL. 

0 LB  WARHEAD*  KK  KILL. 

?0HM  HE  PROJ.*  CLEAN  aircraft*  A KILL. 

?0MM  HE  PROJ.*  aircraft  WITH  WING  TANKS*  A KILL. 

30MM  he  PROJ.*  CLEAN  AIRCRAFT*  A KILL. 

30MM  he  PROJ.*  aircraft  with  wing  tanks*  a kill. 

40MM  HF  PROJ.*  CLEAN  AIRCRAFT*  A KILL. 

40MM  HF  PROJ.*  aircraft  WITH  WING  TANKS*  A KRL. 

5^Mm  HF  PROJ.*  clean  aircraft*  A KILL. 

SIMM  HF  PROJ.*  AIRCRAFT  WITH  WING  TANKS.  A KILL. 

TOTAL  PROJECTED  AREA  OF  A FOREIGN  AIRCRAFT  (FOw  HIT  PKOH.  CaLC.). 
GLAADS  target:  4 SO  M FRONT  K BACK;  ?4  SO  H SIDES*  TOP  K ‘^OTTOP. 


integer  calir.sysin 

DIMENSION  TV(3)  *SVn)  *V(3)  *VR(3)  *RSULT<3*3)  *FR(3)  *«(3)  . 

1S<  7*6.7) *AVR( 10) 


COMPUTE  data  set  NUMBER  AND  READ  DATA  FROM  DISK. 


SYSIN  =10^CALIB 
rewind  SYSIN 

READ  (SYSiN.BOl ) EMAX*  ( ( ( S ( J . K * L ) t J= 1 * 3 ) *K * I * 6 ) * L* I * FM A X ) 
BOl  F0RMAT(I5/(6F10.5) ) 


EMAX  EQUALS  THE  NUMBER  OF  VULNERABLE  COMPONENTS 

K*lJ  FRONT*  Ks2:  REAR*  K=3:  PORT*  K»4:  STaRROAKO*  Kz^:  MOTIOMI 
K*6I  TOP. 

THE  AREAS  IN  THE  ^S*’  MATRIX  ABE  IN  SO  MFTFRS 

CORRESPONDING  VEIOCITIES  FUR  VARIOUS  LEVELS  OF* J aRF  IN  mETFHS/S 


RETURN 


entry  VUI.NARIPOS*  VP.AV) 
real  P0S(6) *AVt7) 
dimension  IN(3) * IM(3) 

p ■ SORT  <POS(n»PnS(l)->POS(?)*PnS(2)^POS(3)*POS(3)  ) 
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M 


FOPTRAN  IV  (5  LEVFL  21 


VULAPn 


OATF  = ^A109 


IS/IS/SS 


OOU 

0015 

0016 
0017 
00l« 
0019 
00?0 
0021 
0022 

0023 

0024 

0025 

0026 
0027 
002R 

0029 

0030 

0031 
003? 


loo33 

0034 

0035 


0036 

0037 
003R 


GP=  50RT  (P05U)^POS(1)^PO«*(2)*POS(2)  ) 

CE  * GP/R 
5E  s P0S{3)/R 
SA  s P0S(2)/r,P 
CA  a P0S(1)/GP 
$V(n  = VP*CE*Ca 
5V(2)  a VP*CE»SA 
SV(3)  = VP*SF 
VR(1)  a SV(1)-P0S(4) 

VR(?)  a SV(2) 

VP(3)  = SV(3) 

PT  a SOPT(  V»n)«VRn)4VP<?)«VP(2)4VR(3)*VR{3)  ) 

DO  14  I = lf3 

13  R(I)  a AHS  (VR(T)/RT) 

AVR(i)  s a6S(  v»{T)  ) 

14  1M(1)  s 5IGN(1.0*VW(I)  )-l.)/2. 

no  10  F a IftMAX 

AV(F)  = 0. 
no  in  1 a 1,1 
c 

C THIS  fOllATtON  I'FPFOMMS  llN^'\P  I MT  f PPL.O  A T T ON  ON  VUl.NFPAHl  f ftOl  AS 

C ANO  S(l.IM(n,F)  M)M  STRlKir^O  VHOClTlrS  oy  V?  ANh  V| 

C PESPFCTI VfcLY.  V?=304.R  M/b  (1000  FT/bFC)l  Vlalb2,4  M/S  (500  F I /Sf r ) , 

C 

ATPa(S(?,IM(I) *F)-S(lt IK (T) tF) )*(AVP(l)-l5?.4)/152.4^S(lf IM(I) tf ) 

ATP  . AMAXI (ATPf 0.) 

ATP  a AMTNl (ATPtS(2f TM(I) ,M ) 

C 

C NOTF  that  linear  I NTFRPl  OLAT ION  IS  USED  TO  COMPUTE  VULNERABLE 

C AREAS.  WF  DO  NOT  ALLOW  an  aRF A SMALLER  THAN  THE  SMALLEST  VAI UF 

C GIVFN  OR  A VALUF  LARGER  THAN  THE  LARGEST  VALUE  GIVEN.  OOn059B0 

C 

30  aV(F)=R(  T)«»ATP^AV(E) 

RETURN 

END 


# 
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SAMPLE  PROBLEM 


The  following  sample  problem  illustrates  the  use  of  the  ISO-PK  pro- 
gram for  the  analysis  of  a hypothetical  air  defense  gun  system  being 
used  against  a hypothetical  target  aircraft. 


Next  page  is  blank. 


> 


1 


INPUT  LISTING 


//FTOlFOOl  DO  SYS0UT=A»DCB=RECFM=UA  OUTPUT  DATA  SET  FOR  FOOTPRINTS 
//DATA  00  * 

3 4 3 T F 

-2000.  4000.  0,  5000.  0.  5000* 

-120,61111-180.05555-231,50000 

250.  450.  050.  1200. 

250.  450.  800. 

SSTOP 

//FT02F001  00  * GUN  DATA  SET 
ISO-PK  SENSITIVITY  ANALYSIS:EXPERIMENTI t 1/2 


FPC/2 

FPC/4 

FPC/5 


1 

12.  .003  .003  1014, 

,60  .60  0.0  30.0  1750.0  1.48 

1500.  2.834 


0.0030 

0.0035  4 

4 

0. 

0.060  0,300 

0.600 

,004  ,0047 

.0058  .0072 

0. 

0.060  0,300 

0.600 

,004  ,0047 

.0058  .0072 

-5, 

5.0 

,0025 

//FTllFOOl  DO  * 
1 

0. 

.0025 

0, 

5,0 

0. 

.0025 

0.  0. 

.0025 

2022110 

.17 

.33 

.33 

.06 

.12 

.12 

2022120 

1.27 

2.54 

2.54 

1.27 

2.54 

2.54 

2022130 

1.11 

2.23 

2,23 

.8 

1.6 

1.6 

2022140 

GADES  rSO-PK  CONT.OUR  PBOGR’Am 
ISO-PK  SENSITIVITY  ANALVSIS:EXPERIMENTI * 1/2 


MARCH 


19  7? 


OUTPUT  LISTING 


caliber  type 

MUZZLE  VELOCITY.  M/S 

STO.  DEV.  OF  MUZZLE  VELOCITY.  M/S 

STO.  DEV.  OF  X-COMP.  OF  RES.  GUN  OISP..  RAO. 

STO.  DEV.  OF  Y-COMP.  OF  RES.  GUN  DISP..  RAO. 

MAXIMUM  ELEVATION.  PAD. 

maximum  elevation  Rate,  rao./sec. 

MAXIMUM  AZIMUTH  RATF*  RAO./SEC. 

AVERAGE  system  REACTION  TIME.  SEC. 

rounds  per  BURST 

MAXIMUM  EFFECTIVE  RANGE.  METERS 

RANGF.  M.  »1SOO.OOOO  TIME  OF  FLIGHT.  SECS 


I 

1014.00 

12.0000 

0.003000 

0.003000 

1.4800 

O.f.0 

0.60 

0.0 

30. 

1750. 

» 2.8340  BALLISTIC  COEFF  TC I ENT  . '<-BAH 


dynamic  gun-pointing 
point  el.  rate.  R/S 
1 0.0 
2 0.060 

3 0.300 

4 0.600 


error  function  - elevation 

STD  OEV  OF  EL.  ERROR.  R 
0.0040 
0.0047 
0.0058 
0.0072 


BIAS.  RAO.s  0.0035 


O' 


dynamic  gun-pointing 

POINT  AZ.  RATE.  R/S 
1 0.0 
2 0.060 

3 0.300 

4 0.600 


ERROR  FUNCTION  - AZIMUTH 
STD  OEV  OF  AZ.  ERROR. 
0.0040 
0.0047 
0.0058 
0.0072 


SENSOR  ERRORS 

MEAN  MEASURENT  ERRORS  XH(I) 


R 


BIAS.  RAO.>  0.0030 


RANGE. M 
AZIMUTH. RAO 
elevation. RAO 
RANGE  RATE. M/S 
AZ.  rate. RAO/S 
EL.  rate. RAO/S 


-5.00000 

0.0 

0.0 

0.0 

0.0 

0.0 


matrix  of  standard  deviations  and  CORRELATION  COEFFICIENTS 
the  units  are  the  Same  as  for  the  means 


RANGE 

AZIMUTH 

ELEV. 

R.  RATE 

AZ.  RATE 

EL.  RATE 

RANGE 

5.00 

0.0 

0.0 

0.0 

0.0 

0.0 

azimuth 

0.0 

0.0025 

0.0 

0.0 

0.0 

0.0 

elevation 

0.0 

0.0 

0.0025 

0.0 

0.0 

0.0 

R.  Rate 

0*0 

0.0 

0.0 

5.000 

0.0 

0.0 

az.  pate 

0.0 

0.0 

0.0 

0.0 

0.0025 

0.0 

EL.  RATE 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0025 

n,32?l 


COWTOUR 

PK 

charactfr 

1 

O.OObO 

F 

2 

0,01U0 

G 

3 

0.0500 

H 

4 

0*1000 

I 

5 

0*2500 

J 

6 

0*5000 

K 

7 

0*7500 

L 

fl 

0*9000 

M 

9 

0*9500 

N 

10 

0*9900 

0 

mask 

character 

EL,  RATE  OR  MAX 

EL 

F 

az.  rate 

A 

INSUF,  TIMF 

T 

OUT  OF  RANOE  R 


5.00  ♦ 


ISO-PK  CONTOUR  MAP 


4.50  • ♦ ♦ 


4,00  * ♦ 


3.50  * ♦ ♦ 


3.00  * 


00 


CROSSING  RNG.f  M 
MULT.  BY  10«#-3 

?.50  * 


?.00  ♦ 


1.50  * 


1.00  * 


0.50  * 


-0.00  

I 

-2.00 


♦ ♦ 


♦ ♦ ♦ ♦ 


♦ ♦ 


RRRRRRRRRRR 

RRRRRR6  6G6RRRRR 

RRR6  GGRRRR 

♦ RRRG  ♦ ♦ ♦ ♦ ♦ G§  RR  ♦ 

RR6  HHHHHHHH  G6  RR 

RRG  HHHHH  HHHHH  G RR 

R6G  HHH  HHH  G P 

RG  HHH  I II II  III  H 

PG^  HHW  mil  ♦ ♦IIII^  HH  ♦ 


SJJJ  ♦ 

JJ 

JJ 

J A 

JJ  KAAAKKKK 


^J 


RG  HH  nil 
RG  H II 
R6G  HH  II 

G H I JJJJJJJ 

♦ RF6  ♦H  !♦ 

F6  H I 

FG  HI 

R F 6 HI 

F 6 HI 

♦-H-f  ♦-d  — ♦KKKIrK K4- 

... I. I*. I.. 


II 


I 


I 

I 

l4 

I 

II 
I 


H 

H 

H 

H 

H 

HH 

H 

HH 


H 

♦ H- 

I .. 


G 

♦ 6 

G 

G 

G 

G 

♦ G 

G 

G 

GG 

6 

-♦6— 


♦ ♦ ♦ ♦ 


RR 

FR 

FR 
F R 
F 

F F 
♦ F 
F 
FF 
F 
F 

-♦p  — 

^ ♦ 


♦ ♦ 


♦ ♦♦♦♦♦♦ 


♦ ♦ 


•1.40 


-0.50 


-0.20  0.40  1.00  1.60  2.20 

DOWN  RANGEf  M.  MULT.  BY  10**-3 


2.80 


...... I . 

3.40 


4.00 


ISO-PK  SENSITIVITY  ANALYSIS:F?tPERIMENTI  1 1/2 
TARGET  VELOCITY*-120.61lM/S  TARGET  ALTITUDE*  250.000m 


♦ 

- 

4 

• 

s.oo 

• _ « 

1 

1 



- - 

ISO-PK  CONTOUR  MAP 





4 

-4- 



— 4- 

1 

1 

4.50 

1 

1 

• 

1 

1 

♦ ♦ 

♦ 

♦ ♦ ♦ 

♦ ♦ ♦ ♦ ♦ 

4 

4 « 

4 

4 

4 

4 

4 

1 

1 

1 

4 

1 

1 

4.00 

1 

1 

» 

1 

1 

♦ ♦ 

♦ 

♦ ♦ ♦ 

♦ ♦ ♦ ♦ 4 

4 

4 4 

4 

4 

4 

4 

4 

1 

1 

4 

1 

1 

3.50 

1 

1 

• 

1 

1 

♦ ♦ 

♦ 

♦ ♦ ♦ 

4 4 4 4 4 

4 

4 4 

4 

4 

4 

4 

4 

1 

1 

4 

1 

1 

3.00 

1 

1 

• 

1 

♦ ♦ 

♦ 

♦ ♦ ♦ 

4 4 4 4 4 

4 

4 4 

4 

4 

4 

4 

4 

1 

1 

4 

1 

CROSSING  RNG.f  M 

1 

1 

1 

1 

MULT. 

BY  10**-3 

1 

1 

• 

1 

1 

2.50 

« 

♦ « 

♦ ♦ ♦ 

4 4 4 4 4 

♦ 

4 ♦ 

4 

4 

4 

4 

4 

1 

4 

1 

( 

VO 

2.00 

1 

1 

1 

« 

1 

♦ ♦ 

♦ 

♦ ♦ ♦ 

4 4 4 4 4 

4 

4 . 

4 

4 

4 

4 

4 

1 

1 

1 

4 

1 

1 

1 

RRRRRRRRRRRPRRRRRR 

1 

1 

1 

PPRRG  GGRRRR 

1 

1.50 

• 

♦ ♦ 

♦ 

♦ ♦PRPG^ 

4 4 4 4 GGIVRR  ♦ 

4 4 

4 

4 

4 

4 

4 

4 

1 

RRG 

GG 

RR 

1 

1 

RRG  HHHHHHHHHHHHHH 

G RR 

1 

1 

RGG  HHHH 

t HHH 

G6  PR 

1 

1 

PPG  HHH 

nil  HH 

G 

R 

1 

1.00 

# 

♦ ♦ 

♦ P6«  HWH  Ullli  mil  ♦ W 

♦ G 

H 4 

4 

4 

4 

4 

4 

4 

1 

RGB  HH  II  II 

II  H 

G 

R 

1 

1 

G6 

HH  II 

II  H 

G 

FP 

1 

1 

RG 

H II 

I H 

G 

F P 

1 

1 

RGG 

HH  II 

JJJJJ  I H 

G 

F 

1 

0.50 

• 

♦ ♦ 

♦ G 

♦ H !♦  J<(J  ♦ !♦  H ♦ 

♦ G 

♦ F P4 

4 

4 

4 

4 

4 

4 

1 

RF6 

H I JJ 

1 J 1 H 

G 

F 

1 

1 

F6 

H II  JJ 

J I H 

G 

F 

1 

1 

F6 

H I JJ 

A J I H 

G 

F 

1 

1 

F6 

H I JJ 

AAAKKKKK  j I H 

6 

F 

1 

-0.00 

F6-- 

KKPf KK f — H-4- 

1 . 

J • 

-2 

.00 

-1.40 

-0. 

BO  -0.20 

0.40  i.oo 

1.60 

2.20 

2.80 

3.40 

O 

O 

. 

DOWN  RANGEf  M.  MULT. 

BY  10** 

-3 

ISO-PK 

SENSITIVITY 

ANALYSIS: 

EXPERIMENTI«1/2 

TARGET 

VELOCITY*' 

-12B 

.611M/S 

TARGET  altitude* 

450.000H 

5.00  ♦ 


ISO-PK  CONTOUR  MAP 


00 

O 


4,50  * 


4,00  ♦ ♦ ♦ ♦ ♦ ♦ 


3.50  * ♦ 


3.00  * 


CROSSING  RNG.f  M 
MULT.  BY  10##-3 

?.50  ♦ 


?.00  ♦ 


1.50  ♦ 


1.00  • 


0.50  * 


•0.00  

I «•. 

-2.00 


♦ ♦♦♦♦♦ 


♦ ♦ 


♦ ♦ ♦ ♦ ♦ 


♦ ♦ 


RRRRPRRRRRRRRR 
♦ RRRR#;  ♦ ♦ GGRRRRR  ♦ 

RRRG  GGGRP 


RPG 

RRG  HHHHHHH 

PRG  HHHHHH  HHHH 

♦RG  ♦ HHM  ♦ ♦ ♦HH  ♦ 

PG  HHH  I III  II  1 1 II  I HH 

RG  HH  III  II  H 


ggrpr 
6 

GG 


RG 

GG 


HH  III 
HH  II 


II 


I 


♦ RG  ♦ HH  ♦II 
GG  H I I 

PG  H 1 1 

G HH  I 

G H I 

♦H  — ♦♦ <^ddddiR 


RR 

P 

G«  P 
G w 
G f 
G 

G 

♦ G 


R 

P 

♦ r 


AA  UUUvJvJU 
JJAAAJJ  J 

♦ 


d-*- 


H 

H 

-M4- 


G P 

G r 

6 F 

0 F 


♦ 4 


♦ ♦ 


-1,40 


-O.RO 


•0.?0  0,40  1.00  1, 

DOWN  RANGEf  M.  MULT,  RY 


. I 

60 

10 


> I . 


I 


I , 


2.20 


P.flO 


3.40 


4.0 


•♦-3 


ISO-PK  SENSITIVITY  ANALYSIS;EXPERIMENTItl/2 
target  VEL0CITY»-I28.611M/S  TARGET  ALTITUDE* 


BOO. DOOM 


ISO-PK  CONTOUR  MaP 

I I 

I 1 

I I 

I I 

4«SO  * ♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦ 

I I 

I I 

I I 

I I 

4«00  * ♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦ 

I I 

I I 

I I 

I I 

3*50  ♦ ♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦ 

I I 

I I 

I 1 

I 1 

3,00  * ♦♦♦♦♦♦*♦♦♦♦♦♦♦♦♦♦♦♦♦ 

I I 

CROSSING  RNG.f  M 1 I 

MULT.  BY  10*«-3  I I 

I I 

2,50  * ♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦ 


2.00  • 


1.50  • 


RRRRRRRRRHR 

RRRRRRGGGGGG6666  RRRRRR 
RRRGG6G  G6GG  RRR 

♦ RRGGG  ♦ ♦ ♦ ♦ G^G  RRR  ♦ 

RRRGG  GG  RR 

RGG  HHHHHHHHHHHH  6 F« 

RRGG  HHHH  HHH  G FRR 


1 

RGfi 

HHH 

HHH 

6 

F 

R 

1 

1.00  • 

♦ ♦ 

♦ *R6G  *HHH  ♦ iiniiimi  ♦ 

♦ H ♦ 

♦ 

G ♦ 

F 1 

R4  ♦ 

4 

♦ ♦ 

♦ ♦ 

1 

RGG  HH 

nil  III 

H 

G 

F 

R 

1 

I 

RGG 

HH 

III 

II 

H 

6 

F 

R 

1 

1 

RGG 

HH  II 

I 

H 

6 

F 

R 

1 

1 

FF 

H II 
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INTRODUCTION 


The  expressions  used  to  define  the  position  of  the  gun  necessary  to 
score  a hit  on  the  target  are  derived  in  this  appendix.  The  equations 
are  based  upon  linear  prediction  methods  and  are,  therefore,  correct  only 

if  the  target  is  flying  at  a constant  speed.  The  expressions  must,  of 

• • • 

necessity,  be  in  terms  of  the  measured  quantities  (a,  e,  r,  a,  e,  r) 

(see  Figure  1,  p.  7)  at  the  present  position  of  the  target;  for  only  then 
can  the  necessary  partial  derivations  be  taken. 

Azimuth  Angle  of  the  Gun 

The  azimuth  position  of  the  gun  (and  target)  is  a rather  artificial 
quantity  since  there  usually  is  no  reference  on  the  gun  from  which  it  can 
be  measured.  However,  it  is  a required  mathematical  necessity.  It  will 
be  shown  that  the  actual  measurement  of  an  azimuth  angle  is  not  required 
in  determining  the  lead  angle  of  the  gun.  The  position  of  the  gun  can 
be  expressed  as: 

Ag  = a + Xa 

where 

Aj,  = azimuth  angle  of  the  gun 

= azimuth  lead  angle  measured  in  the  ground  plane 
Mathematically,  Aq  and  "a"  are  both  measures  from  the  same  arbitrary 
reference. 
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the  lead  angle  is  defined:  Consider  the  projection  of 

the  path  of  the  target  into  the  ground  plane  as  shown  in  Figure  A-1.  The 
XY  plane  defines  the  arbitrary  coordinate  system  in  which  the  azimuth 

angle > a,  is  measured.  The  xy  plane  is  not  arbitrary;  however,  it 
is  set  up  so  chat  the  y-axis  coincides  with  the  reference  position  on 
the  gun  from  which  the  lead  angle  is  measured.  This  axis  also  lines  up 
with  the  measured  position  of  the  target  (T(t)).  The  velocity  vector  v 
to  signify  the  path  of  the  target,  and  the  position  of  the  target  at 
the  time  of  predicted  intercept  (T(t+At))  is  shown  in  Figure  A-1. 

In  the  xy  plane  and  in  terms  of  the  measured  quantities, 

• • 

X « ra  cos(e) 

yo*  r cos(e) 

* • • 

y = r cos(e)  - re  sin(e) 

The  lead  angle  X^  can  be  expressed  as 

■ »n-X 

At  is  time  of  flight  of  the  projectile  to  the  point  of  intercept. 

It  is  more  convenient  to  leave  the  expression  for  the  lead  angle  in  this 
form.  Note  that  all  of  the  terms  in  this  equation  (except  At)  are  func- 
tions of  the  measured  quantities.  At  will  be  treated  later,  and  will  be 
shown  to  be  a function  of  the  measured  quantities. 

The  azimuth  angle  of  the  gun  is  defined  as 

Ac  - a + tan-1 


92 


The  various  partial  derivatives  of  this  expression  with  respect  to  the 


measured  quantities  can  be  obtained  by  application  of  the  chain  rule  for 
partial  differentation.  Note  that  the  lead  angle  calculation  is  indepen- 


dent of  the  azimuth  angle  measurement,  but  that  the  gun  position  can  be 
wrong  if  there  is  a measurement  error  involved  with  the  azimuth  angle  a. 
This  error  can  be  likened  to  a tracking  error  in  the  azimuth  direction. 

Elevation  Angle  of  the  Gun 

The  elevation  angle  of  the  gun  is  shown  in  Figure  A-2.  The  u-axis 
in  this  figure  coincides  with  the  one  shown  in  Figure  A-1.  The  expression 
for  the  elevation  of  the  gun  (Eq)  is  given  by 


Where 


z(t  + At)  = the  altitude  of  the  target  at  time  of  predicted 


intercept , 


R£  = the  range  to  the  target  at  time  of  predicted  intercept 


In  terms  of  the  measured  quantities: 


z(t  + At)  * Zq  + zAt 


Where 


Zq  = r sin(e) 
z = re  cos(e)  + r sin  (e) 


and  (referriag  to  Figure  A-1) 


r|  - (iAt)^  ^0 


- (x2  + y2  ^ ^^2  + 2 rr  At  + r2 


93 


4 


0) 


Figure  A-1.  Ground  Plane  Projection  Of  An  Aircraft  Flying  At  A Constant  Speed. 


+ zAt 


I 


' Figure  A-2.  Elevation  View  Of  The  Target  At  Time  Of  Predicted  Intercept 
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The  terms  defining  Eg  can  be  expressed  in  terms  of  the  measured  quan- 
tities; therefore,  the  necessary  partial  derivatives  can  be  obtained  using 
the  chain  rule.  Note  that  the  actual  measurement  of  the  present  eleva- 


tion angle  of  the  target  is  required  here  to  calculate  the  elevation 
angle  of  the  gun. 

Time  of  Flight 

The  time  of  flight  of  the  projectile  to  the  point  of  predicted  inter- 
cept must  also  be  expressed  in  terms  of  the  six  measured  quantities.  This 
is  required  because  the  partial  derivatives  of  this  quantity  are  necessary 
to  complete  the  total  partial  derivatives  of  the  expressions  for  Aq  and 


E, 


'G* 


This  derivation  is  based  upon  the  two  assumptions: 


1.  Linear  prediction  (i.e.,  the  target  is  not  accelerating) 


2.  The  projectile  obeys  the  following  equation 


R 


■f 


At  * 


Vm  - KRf 


The  derivation  proceeds  as  follows; 


Solving  the  Rf  from  the  equation  for  At 

At 


1 + KAt 


Referring  to  Figure  A-1 


R^  = At^  + 2rr  At  + r^ 


and  equating  to  obtain 


= At^  + 2rr  At  + r^ 


(1  + K At)2 
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Carrying  out  the  Indicated  multiplications  and  combining  coefficients  of 
powers  of  At  yields 

At^Cic^V^)  + (2  rr  + 2v2K)At^  + (r^R^  + 4 rr  K + - V^)  At^ 

+ (2  r^K  + 2 rr)  At  + r^  = 0 

where 

K = ballistic  coefficient  of  the  projectile 
= velocity  of  the  target 
Vm  = muzzle  velocity  of  the  projectile 

The  only  term  in  this  quartic  equation  which  must  be  expressed  in 
terms  of  the  measured  quantities  is  Vj.,  the  velocity  of  the  target.  This 
can  be  expressed  as 

+ y2  + 

= (ra  cos(e))2  -r2e2)  + 

This  completes  the  derivation  of  the  equations  defining  the  position  of 
the  gun  in  terms  of  the  measured  quantities.  The  partial  derivatives  of 
the  expressions  for  and  Eq  can  be  obtained  by  application  of  the  chain 
rule  for  partial  differentation.  Similarly,  the  partial  derivative  of 
At  with  respect  to  the  measured  quantities  can  be  obtained  and  used  in 
the  expressions  for  and  Aq. 

Angular  Rates  of  Guns 

The  angular  rates  of  the  gun  in  the  azimuth  and  elevation  directions 
are  required  so  checks  in  the  ISO-PK  program  can  be  made  against  maximum 
rates  allowed.  These  rates  can  be  obtained  from  the  expression  derived 
earlier  in  this  appendix. 
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The  partial  differential  equations  necessary  to  obtain  F on  page  13 
are  contained  In  the  section  titled  FORTRAN  Source  Program  Listing  on 
pages  56,  57,  and  58  of  the  listing  in  subroutine  GFUNC.  The  deriva- 

tion of  these  partlals  is  straight  forward  but  very  tedious.  It  was  felt 
that  inclusion  of  the  derivation  of  these  partlals  would  be  of  little 
interest  to  readers  or  users  of  this  report.  Therefore,  only  the  final 
form  of  the  partlals  in  FORTRAN  notation  is  included.  The  terms 
are  given  FORTRAN  names  which  make  a mnemonic  interpretation  of  the 

quantity  easy.  For  example,  PAXD  stands  for  the  partial  of  azimuth 

3A 

with  respect  to  2^-^ot  or  — . 

3x 


i 
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DRSAR-RDG 
DRSAR-ASV 

Rock  Island,  IL  61201 
Commander 

US  Army  Test  and  Evaluation  Command 
ATTN:  DRSTE-SY 

Aberdeen  Proving  Ground,  MD  21005 
Commander 

US  Army  Electronics  Command 

ATTN:  DRSEL-SA 

Fort  Monmouth,  NJ  07703 

Commander 

US  Army  Missile  Command 

ATTN:  DRSMI-CS 

Redstone  Arsenal,  AL  35809 

Commander 

US  Army  Tank  Automotive  Command 
ATTN;  DRSTA-V 
Warren,  MI  48090 

Commander 

HQ,  US  Army  Aviation  Command 
ATTN:  DRSAV-D 

P.O.  Box  209,  Main  Office 
St.  Louis,  MO  64502 
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1 Commander 

US  Army  Troop  Support  Command 

ATTN:  DRSTS-G 

St.  Louis,  MO  63120 

1 Commander 

US  Army  Materiel  Systems  Analysis  Activity 
ATTN:  DRXSY-AD 

Aberdeen  Proving  Ground,  MD  21005 
Commander 

Rock  Island  Arsenal 
6 ATTN:  SARRI-LR 

2 SARRI-LPL 

Rock  Island,  IL  61201 

3 Commandant 

US  Army  Air  Defense  School 

Fort  Bliss,  TX  79916 

1 Project  Manager  for  Army  Radar-Directed  Gun  ' 

Air  Defense  System 
US  Army  Armamant  Command 
ATTN:  DRCPM-ARGADS 

Rock  Island,  IL  61201 

1 Commander 

Frankford  Arsenal 
ATTN:  SARFA-FCW-N4200 

Philadelphia,  PA  19137 

1 Commander 

US  Army  Aviation  Center 
Fort  Rucker,  AL  36362 

^ Office  of  Project  Manager  for  XMl  Tank  Systems 

ATTN:  DRCPM-GCM-SI 

Warren,  MI  48090 

12  Defense  Documentation  Center 

Cameron  Station 
Alexandria,  VA  22314 
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